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This  thesis  is  an  application  Df  parameter  sensitivity 
analysis  to  aerodynamic  parameters  of  a  Bank-to-Turn 
missile.  In  the  development  a  brief  review  of  trajectory 
sensitivity  theory  is  presented.  A  linear  analysis  is 
performed  using  an  Uncoupled  Pitch  Channel  Autopilot  and  a 
Coupled  Roll-Yaw  channel  Autopilot  of  the  missile  taken  as 
model.  Finally, a  nonlinear  analysis  is  given  for  the  system. 
Comparisons  between  the  linear  and  nonlinear  cases  are 
outlined. 
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roll  attitude  command  from  guidance  computer,  zero  degrees 
in  -z  direction  and  90  degrees  in  y^direction 

roll  attitude,  zero  degrees  in  -redirection  and  90  degrees 
in  y^  direction 

roll  attitude  error. 

Elevation  Euler  Angler,  second  rotation,y^(qco3  p  -  rsin  <|)  )  dt 

Azimuth  Euler  Angle,  first  rotation 
about  y^  ,  J" (gsin  (p  +  rcos  p  )  dt 

pitch  control  incidence  (positive  tail  incidence  produces 
negative  pitiching  moment) 

commanded  pitch  control  incidence, 

yaw  control  incidence  (positive  tail  incidence  produces 
negative  yawing  moments) 

commanded  yaw  control  incidence, 

roll  control  incidence  (positive  tail  incidence  produces 
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positive  rolling  mcment) 
commanded  roll  control  incidence 
constant  or  equilibrium  angle-of-a ttack 
angle-of-attack 
angle-of-attack  rate 

modified  form  of  estimated  angle-of-attack  for  autopilot 
coordination  command 

angle  of  sideslip 

sideslip  angular  rate 
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I-  parameters  sensitivity  analysis  op  a  BANK-TO-raRN  MISSILE 

The  determination  of  changes  in  system  performance  due 
to  changes  in  parameters  is  of  great  importance  in  engi¬ 
neering  analysis  and  design.  Sensitivity  guest  ions  arise 
when  model  uncertainty  is  present,  or  when  a  range  of  oper¬ 
ating  conditions  is  contemplated. 

The  guestions  of  parameters  sensitivity  particularly  arise 
in  the  fields  of  engineering  where  models  are  idealized, 
inexactly  identified,  or  the  systems  themselves  are  subject 
to  unpredictable  changes  with  tiae  due  to  en viromental, 
material  property  or  operational  influences  so  that  there  is 
always  a  discrepancy  between  the  physical  reality  and  the 
mathematical  model. 

Sensitivity  analysis  provides  the  engineer  with  methods  for 
investigating  or  minimizing  such  parameter  deviations. 

In  general,  the  diagram  of  a  system  can  be  represented 
by  a  single  block  as  given  on  Fig.  1.1. 


SYSTEM 

Figure  1.1  General  Representation  of  a  Dynamic  System. 

From  a  mathematical  point  of  view,  what  we  call  a 
system  is  the  explicity  or  implicity  given  relationship 
between  the  input  signal  u  and  output  signal  y.  In  general, 
u  and  y  can  be  vectors.  The  character  of  this  relationship 
is  commonly  called  the  structure  of  the  system. 
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For  example,  the  structure  of  the  system  may  be 
characterized  by  the  order  of  a  differential  or  difference 
equation,  linearity  or  nonlinearity,  the  order  of  the  numer¬ 
ator  and  denominator  of  a  rational  transfer  function  and  the 
quantitative  properties  of  the  system  parameters. 

Typical  parameters  are  initial  conditions,  time- 
invariant  or  time-variant  coefficients,  natural  frequencies 
and  samplinq  periods. 

The  chanqe  of  the  state  or  the  change  of  the  output 
variable  with  time,  can  be  caused  by;  (1)  the  influence  of 
input  signals,  (2)  the  change  of  parameters.  These  quanti¬ 
ties  are  shown  in  Fig.  1.2. 


Figure  1.2  Quantities  Affecting  the  Dynamics  of  a  System* 

This  work  addresses  the  application  of  the  sensitivity 
analysis  to  the  variation  of  aerodynamic  parameters  of  a 
Bank-to-Turn  missile. 

PARAHETER  SENSITIVITY  is  the  affect  of  parameter  changes  on 
the  dynamics  of  a  system,  say,  the  time  response,  the  state, 
or  any  other  quantity  characterizing  the  system  dynamics. 

In  order  to  have  a  realistic  model  to  work  with  a  NASA 
Contractor  Report  [Ref.  2]  was  adopted  as  a  reference  for 
application  of  the  sensitivity  analysis  on  the  variation  of 
aerodynamic  parameters  of  a  Bank-to-Turn  missile.  Appendix 
A,  B  and  C  give  the  detailed  aerodynamic  assumption  of  the 
system. 
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One  brief  explanation  of  the  sensitivity  theory  and 
analysis  is  given  as  a  guideline  to  better  understand  the 
subsequent  work. 

An  analysis  of  the  sensitivity  of  the  aerodynamic  parame¬ 
ters  applied  to  the  case  of  the  linear  uncoupled  pitch  and 
coupled  roll-yaw  autopilot  is  performed. 

The  pitch  and  roll- yaw  autopilots  are  coupled  as  given 
in  Fig.C.1  in  the  appendix  C.  A  nonlinear  parameter  sensi¬ 
tivity  analysis  is  performed. 

Conclusions  about  the  influence  of  parameters  of  concern 
are  given  and  a  comparison  between  the  linear  and  nonlinear 
case  is  done.  Comments  and  recommendations  for  the  future 
are  outlined  in  order  to  delineate  the  continuity  of  the 
present  work. 
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II.  SENSITIVITY  THEORY 


A.  INTRODOCTIOH 

The  basis  of  all  sensitivity  considerations  in  the  case 
of  time-invariant  parameter  variations  is  the  so-called 
sensitivity  functions.  Dynamic  systems  can  be  characterized 
in  several  ways: in  the  time  domain,  in  the  frequency  domain, 
or  in  terms  of  a  perfomance  index.  There  is  evidently  an 
adequate  number  of  ways  to  define  the  sensitivity  functions 
of  a  dynamic  system.  The  definition  that  is  actually  used 
depends  on  the  form  of  the  mathematical  model  as  well  as  on 
the  purpose  under  consideration. 

The  sensitivity  functions  can  be  classified  into  the 
following  three  categories: 

(1)  Sensitivity  functions  in  the  time  domain 

(2)  Sensitivity  functions  in  the  frequency  domain,  and 

(3)  Per fomance- index  sensitivity 

In  this  chapter  we  will  outline  just  the  sensitivity  func¬ 
tions  in  the  time-domain  of  the  continuous  systems. 
Aplications  of  this  analysis  follow  throughout  this  work. 

The  reader  can  find  detailed  informations  about  the 
others  two  categories  of  sensitivity  functions  in  the 
[Ref.  1], 

B.  TRAJECTORY  SENSITIVITY  FUNCTION  OF  CONTINUOUS  SYSTEMS 

A  continuous,  possibly  nonlinear  system  of  nth  order 
can,  in  general,  be  described  in  the  state  space  by  a  vector 
differential  equation  as  seen  in  Eqa.2.1. 
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X  =  f(X,t,U,'<o) 


(2.1) 


,  X(t^)  =  XO 

Here  X  is  an  n  x  1  state  vector,  f  an  n  x  1  vector  func¬ 
tion,  "u  an  input  vector,  a  nominal  r  x  1  parameter 
vector,  and  X  is  the  n  x  1  initial  condition  vector  or 
initial  state.  Eqn.2.1  is  called  the  NOMINAL  STATE 
EQUATION. 

Assuming  that  the  parameter  vector  deviates  from  the 
nominal  value  by  we  have  the  Eqn.2.  2. 


X 


f(X,t,U,  OC),  X(to) 


XO 


(2.2) 


with  the  initial  conditions  XO  unchanged.  This  equation 
is  called  the  ACTUAL  STATE  EQUATION. 

Now  it  is  assumed  that  Eqn.2. 2  has  a  unique  solution,  X= 
X(t,^)  for  all  admissible  initial  conditions  and  parameter 
values. 

X  is  of  course  a  function  of  u  ,  XO  and  t  as  well. 
However  ,  this  dependence  is  not  needed  for  the  following 
considerations  and  will,  therefore,  be  dropped  for  ease  of 
notation.  Furthermore,  the  solution  X  is  assumed  to  be  a 
bounded  continuous  function  in  t  and  ^  . 

If  the  parameter  takes  on  its  NOMINAL  value  <x-o  »  the 
nominal  solution  X^  =  X(t,o^)  is  obtained.  If,  on  the 
other  hand  the  ACTUAL  solution  is  given  by  X  ”=  X(t,  «<  )  , 
then  the  parameter- induced  change  of  the  state  vector  is 
given  by 


JS  X(t,<K)  =  X(t,  X)  -  X(t,  CK^) 


(2.3) 
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k  first-order  approximation  of  ^  X  can  be  obtained  by 
using  Taylor  expansion  in  the  form  of  Eqn.2.4. 


A  X  (t ,  o<  ) 


2 

Jr1 


(2.4) 


This  equation  can  be  viewed  as  a  definition  of  the 
parameter-induced  trajectory  deviation. 

The  subscript  oc©  shall  indicate  that  the  partial  deriv¬ 
ative  expressed  by  ”7)  is  taken  at  nominal  parameter  values. 

Let  the  state  X  of  a  continuous  system  be  a  continuous 
function  of  a  time- invariant  parameter  vector  o<  = 
{ o<,  ,  c<2  •  Then  the  partial  derivative  can  be 

defined  as 

-  - 

A;(t,^«)  =  VS”*”” 

x*"  1  f  2^  •  •  •  /  n 

j— 1 »  2#  .  .  .  /  r 

fiqn.2-5  is  called  the  trajectory  sensitivity  vector  with 
respect  to  the  jth  parameter. 

Note  that  the  trajectory  sensitivity  vector  is  of  the 
same  dimension  as  the  state  vector,  namely,  n.  Its  compo¬ 
nents  are  the  TRAJECTORY  SENSITIVITY  FONCTIONS  as 


T  ..  .  A  (t,«.) 


(2.6) 


Eqn.2.6  is  the  partial  derivative  of  the  ith  state  vari¬ 
able  in  relation  to  the  jth  parameter  as 
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(2.7) 


.X 


^0 


—  .j 

D  <J  OCo 


Hence,  all  n  x  r  trajectory  sensitivity  functions  form 
the  trajectory  sensitivity  matrix  as  given  in  Eqn.2.8  or 

2.9. 


X 


^  6< 


r  "t)  "1)  xi 

-  »  •  •  I— 


A 


^Xn  ^  ^  Xn 


(2.8) 


(2.9) 


The  columns  of  A  are  the  trajectory  sensitivity 
vectors  Aj-  .  Here  X  is  the  Jacobian  matrix  of  the  state 
vector  with  respect  to  the  parameter  vector  «.  ,  taken  at 

the  nominal  parameter  values. 

With  these  definitions  the  parameter-induced  change  of 
the  trajectory  can  be  taken  as 


ir  __ 

Zilct,^)  =  A  (t,  <x^)  =  X  (2.10) 


Where  o<  =  ,  which  is  the  iCTOiL  parameter 

vector  of  the  system. 
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C.  TBAJECTORY  SENSITIVITY  EQUATIONS  OF  CONTINOOOS  SYSTEMS 

Lets  consider  the  general  continuous  system  described  as 
previously  by  the  state  equation  (Eqn.2.1). 


X  =  f  (X,t,U,  X),  X(t^)  =  xo 


(2.  11) 


Where  X  denotes  the  n- dimensional  state  vector,  ^  the 
r-dimensional  parameter  vector,  f  an  n-dimensional  vector 
function,  and  u  the  input  vector  independent  of  c><.  .  it  is 

assumed  that  the  continuity  conditions  are  fulfilled  and 


that  o<.  is  time-invariant 


Considering  c<  -parameters  and  taking  the  partial 
derivative  of  X  (Eqn.2.11)  with  respect  to  «<j'  ,  one  obtains, 
by  the  application  of  the  chain  rule 


(2.  12) 


The  derivative  of  the  initial  conditions  vector  X®  with 
respect  to  c<j‘  is  zero,  since  XO  does  not  depend  on  X  . 

If  X  is  r-dimensional,  there  are  r  equations  of  the 
form  of  Eqn.2.1 2.  If  we  now  interchange  the  sequence  of 
taking  the  derivative  wirh  respect  to  time  t  and®*<’j’,  and 
then  let  X  approach  to  ^  one  obtains  the  following 

equation 


(2.  13) 
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Here  is  the  trajectory  sensitivity  vector 

-J  Q 

with  respect  to  the  jth  parameter. 

The  Eqn.2.  13  is  called  the  sensitivity  equation  in  the 
state  space  or  the  TRAJECTORY  SENSITIVITY  EQOATION. 

The  above  shows  that  for  c<  -parameters  all  initial 
conditions  of  the  trajectory  sensitivity  equations  are  equal 
to  zero. 

Now  consider  the  output  vector  equation,  as  given  by 


Y  =  g(X,t,u,<X.) 


(2.  14) 


In  a  procedure  similar  to  that  above  one  obtains  the 
algebraic  sensitivity  equation  as  seen  in  Eqn.2. 15. 


I- 


i;. . 


-• 


(2.15) 


the  output  sensitivity 


vector  = 
A.- 


_  __  which  relates 

0  y 

to  the  trajectory  sensitivity  vector  ^j'  .  This 
equation'^is  called  the  VECTOR  OUTPUT  SENSITIVITY  EQOATION. 

Using  the  trajectory  sensitivity  matrix  A  and  the 
output  sensitivity  matrix  C*  ,  the  above  result  can  be 
rewritten  in  the  following  general  form 


i:  T 

A  =  — r. — U  A  + 


1)4  T 

1_  .X->  = 

Ooi.  o<o 


(2. 16) 


cr  =  X  ♦ 


^  oe.  ‘Xo 


(2.  17) 
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These  equations  are  called  the  STATE  SENSITIVITY 
EQUATIONS  of  the  system.  It  is  seen  that  these  equations 
are  linear  whether  the  original  system  is  linear  or  nonlin¬ 
ear. 

If,  in  particular,  the  original  system  is  Linear,  the 
state  equations,  take  the  form 


X  =  AX  +  Bu  ,X(to)  =XO 


(2. 18) 


y  =  c  X  +  D  u 


(2.19) 


Where,  in  general, 

r=T(o<.)  ,'b  =  ¥(o<.)  ,c  =  c(o<)  ,d  =  ¥('^), 

"x  =  ?(t,c<  )  ,  amd  7  =  7(t,6<  )  . 

Note,  however,  that  u  is  not  a  function  of  «<.  if  u  is 
defined  as  an  external  input  of  the  system. 

Now  taking  the  partial  derivatives  with  respect  to  , 

reversing  the  order  of  differentiations  with  respect  to 
time  and  ,  and  letting  ^  approach  to  ,  the 

TRAJECTORY  SENSITIVITY  EQUATIONS  are  obtained 


Xj  .  s,  X\  *  ^ - 1  xo  . 


I-  u(t)  ,  ACl:o)  =  0  (2.20) 


Where  A^  =  A(o<.o)  »  =  X(t,  -^To)  ,  and  j=  1, 


2,...,  r. 


The  initial  condition  vector  X  -Lt:©),  is  again  zero  since 

vj 


X(to)  does  not  depend  on  c<  . 
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By  similar  procedure  the  VECTDR  SENSITIVITI  EQUATION 
becomes 

~  --  'D  C  -  'I)  ^  - 

- LXo  ♦  - l_u(t)  (2.21) 

®  lio<^  oc^ 

Where  =  C  (  ck^  )  and  j=1,2,...,  r. 

In  the  case  of  linear  systems,  the  vector  sensitivity 
equations  have  the  same  A  matrix  as  the  nominal  state  equa¬ 
tions  and  hence  the  same  characteristic  matrix  si  -  A. 

They  differ  from  the  nominal  original  state  equation  only 
in  the  driving  function  and  the  initial  conditions.  The 
latter  are  all  zero.  The  driving  functions  can  be  obtained 
by  solving  the  nominal  state  equations. 

Here,  in  the  linear  case,  the  sensitivity  equations  can 
be  joined  to  the  system  equations  ,  forming  the  so-called 
COHBINED  SYSTEM  as  given  in  Eqns.2.22  and  2.23. 


• 

> 

a 

0 

_ 1 

- 1 

CPJ 

0 

• 

"I>  A  1  ^ 

+ 

■^8 

.  X  . 

«<o-‘ 

u.c-t) 


(2. 22) 


■  Co  O 

Xo 

Do 

* 

— 

1 

+ 

'Do 

‘"'’J 

(2.23) 


The  simultaneous  solution  of  these  differential  equa¬ 
tions  using  the  same  standard  procedure  for  each  of  the 
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above  matrix  equations  yields  the  output  vector,  the  state 
vector,  and  the  corresponding  output  and  trajectory  sensi¬ 
tivity  vectors  ,  and  Aj  •  If  there  are  r  parameter 

variations,  r  systems  of  equations  of  the  above  type  have  to 
be  solved.  Since  the  homogeneous  part  of  the  original 
differential  equation  is  identical  with  the  homogeneous  part 
of  the  sensitivity  equations  with  respect  to  all  parameters. 

A  graphical  interpretation  of  Eqns.2.22  and  2.23  is 
given  in  Fig.2.1. 

Since  the  driving  function  of  the  sensitivity  model 
contains  the  nominal  state,  the  measuring  circuit  for  the 
trajectory  sensitivity  functions  consists  of  a  connection  of 
the  nominal  original  system  and  the  sensitivity  model  as 
illustrated  by  Fig.2.1. 

If  the  actual  input  u  is  applied  to  this  structure  the 
trajectory  and  output  sensitivity  vectors  Aj  and  <3^-  can  be 
measured  at  the  points  1  and  2  respectively  in  the  Fig.2.1. 

In  order  to  measure  all  r  sensitivity  vectors  simultane¬ 
ously,  r  sensitivity  models  are  required. 

D.  STROCTURAL  METHOD 

The  basis  for  the  measurement  of  the  trajectory  sensi¬ 
tivity  functions  is  the  STRUCTURAL  INTERPHETATION  Of  the 
trajectory  sensitivity  equation. 

The  physical  system  represented  by  the  trajectory 
sensitivity  equations  is  called  the  trajectory  sensitivity 
model  or  the  sensitivity  model  in  the  stare  space. 

Regardless  of  the  nature  of  the  original  system  the 
sensitivity  model  in  the  state  space  is  always  linear.  The 
graphical  illustration  of  sensitivity  equation  can  be  given 
by  Fig. 2. 2.  The  system  described  by  the  sensitivity  equation 
is  refered  to  as  the  sensitivity  model  of  the  original 
system.  For  each  output,  a  system  has  as  many  sensitivity 
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models  as  parameters  of  interest.  Both  the  nominal  original 
system  and  the  corresponding  sensitivity  models  form  the 
COMBINED  SYSTEM.  The  sensitivity  model  is  always  linear.  If 
Y(t,^)  is  the  output  of  a  syste^  the  corresponding  sensi¬ 
tivity  functions  (t,c<p)  =  outputs  of 

the  corresponding  sensitivity  models.  Consequently,  in  order 
to  measure  the  output  sensitivity  functions  simultane¬ 
ously, the  nominal  original  system  and  the  sensitivity  models 
have  to  be  measured  at  the  outputs  of  the  sensitivity 
models. 

The  double  frame  used  for  the  original  system  in  Fig. 2.2 
is  to  indicate  that,  in  general,  the  system  may  be  nonlin¬ 
ear,  whereas  the  sensitivity  models  are  always  linear.  In 
the  nonlinear  systems  the  Eqn. 2.12  applies  and  different 
programs  have  to  be  set  up  for  the  original  and  the  sensi¬ 
tivity  equations. 

For  some  applications,  the  sensitivity  functions  with 
respect  to  all  parameters  are  required  simultaneously.  If 
there  are  r  parameters  of  interest,  r  sensitivity  models 
would  be  needed.  This  implies  a  rather  extensive  computer 
time . 

A  method  of  determining  all  output  sensitivity  functions 
of  a  system  simultaneously  by  a  single  sensitivity  model  is 
available  ,  which  is  also  called  the  method  of  SENSITIVITY 
POINTS.  This  method  has  application  just  for  linear  systems. 
Detailed  explanation  about  the  application  of  the  sensi¬ 
tivity  points  theory  can  be  found  in  the  [Ref.  1]. 

E.  OOTPOT  SENSITIVITY  PONCTION  OP  CONTINOOOS  SYSTEMS 

Consider  the  input-output  behavior  of  a  continuous, 
possibly  nonlinear,  single-variable  system  described  by  an 
ordinary  differential  equation  of  the  type 
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•  »y»t,  o<^  }  —  0 


(2. 24) 


f{  y 


in)  (n-<; 


with  the  initial  condition  y(tj,)  =  y.  o,  i=0,  1,...n-1.  y 
denotes  the  output  signal,  t  the  time,  and  cx  a  single 
time -invariant  or  slowly  varying  -parameter  that  has 
nominal  value  • 

In  general,  f  is  a  function  of  the  input  u  as  well. 
However,  if  u  is  an  external  input  which  does  not  depend 
on  o<^  ,  the  dependence  of  f  on  u  is  not  relevant  in  further 

considerations , 

Let  us  suppose  that  the  above  NDMINAL  differential  equa¬ 
tion  has  the  unique  solution  given  below 


lo  =  Y(t,  ®<o) 


(2.  25) 


Which  one  shall  call  the  NOMINAL  SOLOTION. 

Let  us  now  assume  that  the  parameter  changes 
from  to  <X.=  ,  where  ^  is  time-invariant  or 

slowly  varying  with  time.  <X.  is  called  the  actual  parameter 
value.  The  corresponding  ACTUAL  DIFFERENTIAL  EQUATION  can, 
then,  be  written  as 

in) 

f  —  {iJ,y  ,...,y,t,oc)—  0  (2. 26 ) 


Note  that  by  this  change  of  o<o  into  oc.  the 
conditions  remain  unchanged,  namely  y(t^)  =  y.®« 


initial 
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The  corresponding  solution  is  given  as 


y  =  y  (t,  o< ) 


(2.  27) 


Which  we  shall  call  the  ACTOAL  (or  PERTURBED)  SOLUTION. 

It  is  assumed  that  y(t,  <x)  is  of  the  same  type  as 
yCt,  «<o  )  »  and  y(t,  o<)  deviates  infinitesimally  from  y(t,o<o) 
if  'X  deviates  infinitesimally  from  •  The  conditions 

for  fulfilling  this  requirement  are  given  in  the  mathemat¬ 
ical  literature.  For  our  purpose,  it  is  sufficient  to  know 
that  y  is  continuous  in  if  f  is  continuous  in  y  which  is 

true  for  all  continuous  systems  and  "Xq  ^  0. 

With  the  above  assumptions,  the  actual  solution  y (t, 

+  ^ae.)  can  be  expanded  into  a  Taylor  series  around  , 

yielding  the  Eqn.  2.28. 


--fl  +. 


(2.  28) 


y(t,<x)  =  y(t,Qg  +  - 


If  Xq  ,  the  Taylor  series  can  be  truncated  at  the 

linear  term. 

This  gives  the  Eqn. 2.29. 


y(t,o^)  =  y(t,(Xo)  +  - - 1  ^ 

()  X  Xo 


(2.29) 


For  finite  values  of  ^  X  ,  this  expression  can  be 
considered  a  first-order  approximation  of  y(t,o<.). 

The  actual  output  can  be  written  as 
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y(t,<x)  =  y(t,c<o)  ♦  (t,o<o)  a  bc 


(2,  30) 


Where  y(t,c*to)  is  the  nominal  output  and  (T  ( 
the  parameter-induced  output  error. 

With  (T  (t,  o<.o)  defined  |  .  Ttie 

induced  output  error  is  in  these  terms  defined  as 


y  (t,^  ) 


A 


^  (tr«o) 


A  OC 


parameter 

(2.  31) 
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Figure  2.1  Sraphical  Interpretatiaa  of  Eqns.2.22  and  2.23. 
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Figure  2.2  Strutural  Diagram  of  the  Combined  System. 
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III.  APPLICAriON  OF  SENSITIVITY  TO  UNBAR  SYSTE MS 

A.  IHTRODOCTION 

The  sensitivity  theory  described  previously  will  be 
applied  to  the  case  of  an  uncoupled  pitch  autopilot  and  to 
a  roll-'yaw  coupled  autopilot  presented  in  the  appendix  C. 
For  each  case,  the  linear  equations  of  the  nominal  system 
are  presented  in  state  variable  form  and  the  correspondent 
trajectory  sensitivity  equations  are,  derived.  A  sensitivity 
analysis  of  the  systems  is  performed. 

B.  ONCOOPLED  PITCH  AOTOPILOT  ANALYSIS 

Linear  Equations  of  the  Nominal  System 

From  the  block  diaqram  of  the  Fig.B.1  in  the 
appendix  B  one  can  obtain  the  following  nominal  state  equa¬ 
tions  of  the  uncoupled  pitch  autopilot. 


XI 


=  C2  A3  X2  +  C2  A4  X3 


(3.1) 


X2  =  XI  -  K  Cl  A1  X2  -  K  Cl  A2  X3 


X3  =  -  C3  X3  +  C3  Conv  X6 


X4  =  -  Cl  C4  A1  X2  -  Cl  C4  A2  X3  -C4  X4 


(3.2) 


(3.3) 


(3.4) 
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X5  =  C7  X4  -  C5  X5  -  C6  NZC 


(3.5) 


X6  =  C9/Conv  +  C2  C8/conv  A3  X2 


(3.6) 


♦  C2  C8/Conv  A4  X3  -  C7  C8  X4  ♦  (C5  C8  -  C9)  X5 
+  C6  C8  HZC 


For  the  purpose  of  using  the  linear  method  these 
state  equations  must  be  presented  in  matrix  form  as  given  by 


o 

o 

o 

•o 

• 

o- 

1 

-K,C,A, 

-KC^A^ 

o 

o 

o 

o 

^3 

o 

o 

o 

o 

CjCOhv 

^3 

+ 

o 

o 

“  ^ 

-^4 

o 

o 

o 

o 

o 

o 

o 

# 

i 

C9 

CrCg  As 

CfcCe  A*q 

-^7^=8 

^6 

6 

1 

conv 

Qonv 

conv 

Nzc 


(3.7) 


The  correspondence  of  the  state  vectors  is: 

XI  =  c|  ,  12  =  ^  ,  13  =  ,  X4  =  X  /  X5  =  Y  ,  X6  =  S 

Definition  of  the  constants  Cl  through  C9  are  given 
in  appendix  B.  The  parameters  of  interest  for  this  system 
are  given  by  A1,  A2,  A3,  and  A4  which  are  : 


A1  =  C.,  ,  A2  =  ,  A3  =  C 

This  nomenclature  is 
sensitivity  theory  and  to 
carrying  symbols  and  constants 


,  A4  = 


used  here  to  easily  apply  the 
avoid  the  inconvenience  of 
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I 


I 


I 


m 


% 


2.  Sensitivity  Equations 

To  apply  the  procedure  developed  in  chapter  2, 
consider  the  trajectory  sensitivity  equation  (Eqa.2.13). 


=  A 


- 

From  Bqn.3.7  one  can  see  that 


25 


J  0&4.;  Zi 


,  X  {%>)zo  (3.8) 


V 


o 

C^Aj 

0 

0 

1 

-H,C,A, 

-KC.A^ 

0 

0 

0 

o 

o 

0 

0 

CjConv 

o 

-CiC^A, 

-^4 

0 

0 

o 

o 

o 

Cf 

0 

conv 

^z<^e  Aj 
conv 

C1C9  A. 
conv 

%  0 

(3.9) 


b-  [  O  O 


-  C 


C 


Ce  Cs  3 


(3. 10) 


Hi  A  D  ft 

The  partial  derivatives  -^-^-.-|^and  -l^^are 

evaluated  considering  the  parameters  of  interest  A1,  A2,  A3 
and  AU,  that  are  respectively  the  aerodynamic  coefficients 
present  in  the  pitch  channel. 
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Applying  the  partial  derivatives  with  respect  to  the 
parameters  we  have  the  following  matrices 


HA 


o 

o 

Ca 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O  ' 

o 

6 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

Ca  C0 

CotxV 

o 

o 

o 

(3.11) 


o 

<^2 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

■JA  . 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

C^CQ 

conv 

o 

o 

o 

o 

(3.12) 


1)  Ajj 


o 

o 

o 

o 

o 

o 

o 

o 

-KC, 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

(3. 13) 


"b  A 
7)  A/ 


o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

<? 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

(3.  14) 
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because  B  is  independent  of  the  parameters  of  interest. 
In  terms  of  components  one  can  see 


(3.  15) 


i=1,2,...,6  j=1,2,3  and  4 

Here,  four  models  are  necessary  to  study  the  effect 
of  parameter  variations.  These  models  are  shovrm  in  Fig.3.1. 

For  instance,  when  one  apply  the  theory  for  one 
parameter  of  interest,  say  A1,  the  sensitivity  equations  can 
be  obtained  from 


\  • 

•  •y. 

r 

^1 

VH 

11 

^2 

^3 

— 

o> 

• 

# 

# 

• 

+ 

^3 

A" 

1 

(3.  16) 


Similar  procedure  can 
ters  A2,  A3,  A4. 


be  done 


for  the  other  parame 


Pigure  3.1  Nominal  and  Sensitivity  Models.  / 
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C.  SENSITIVITY  ANALYSIS 


A  computer  program  using  [Ref.  4]  was  writtsn  in  order 
to  simulate  simultaneously  the  nominal  and  the  sensitivity 
models  as  shown  in  the  appendix  D.  A  step  was  applied  as 
input  of  the  nominal  system. 

The  number  of  equations  solved  are  6  for  the  nominal 
model  and  24  for  the  sensitivity  model.  Each  parameter  was 
varied  simultaneously  1056  from  all  the  nominal  values. 


D.  ANALYSIS  OP  -PARAMETER  VARIATIONS 


The  results  are  plotted  in  Fig. 3.2  to  3.6  and  Table  I 
Each  state  variable  output  is  plotted  separately  with  the 
correspondent  four  output  sensitivity  functions.  By  means  of 
these  plots, the  following  observations  can  be  made: 

The  plot  of  /I,;  =  •  j=1»2,3,  and  4  in  fig.  3.2 

indicates  that  a  parameter  change  /a  Ajg>  ( j=  1 ,2,  3 ,  and  4) 
primarily  affect  the  rise  time  and  overshoot  of  X1(t)  since 
>i11,  /112,  /113,  A14,  are  largest  at  the  time  where  these 


effects  in  XI  occur.  ^ 

The  plot  of  A-2^-  j=1»2,3,  and  4  in  Fig.  3. 3 

indicates  that  X2(t)  is  little  affected  due  to  ^  Aj^  and 
A  A20  and  strongly  affected  in  the  rise  time  due  to  param¬ 
eter  variations  ^Ag^  and  ^  A^c><  ^^e  overshoot  and  steady 

state  are  little  affected  due  to  parameter  variations 


The  plot  of  >l2,  =  »  j=1#2,3,  and  4  in  Fig.  3.4 


/hAj  (j=1,2,3,  and  4) 

''il  'I)  Aj  o  ' 
indicates  that  X3(t)  is  little  affected  in  the  rise  time  due 

to  parameter  variations  A  (j=1,2,3,and  4).  The  overshoot 
and  steady  state  are  little  affected  due  to  parameter  varia¬ 
tions  ^A|o  and  AKzo  and  strongly  affected  due  to  A  Ago  and 


AAj-0. 
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The  plot  of  j  =  j=1»2,3,  and  4  in  Fig.  3.5 

indicates  that  X4(t)  is  little  affected  in  the  rise  time  and 

overshoot  due  to  parameter  variations  ^  Ato 

strongly  affected  due  to  ^  and 

The  plot  of  Atr‘,=  r  j=1,2,3,  and  4  in  Fig.3.6 

5J  A_i  Aj  o 

indicates  that  X5(t)  is  little  affected  in  the  rise  time  due 
to  parameter  variations  A  A^q  (j=1»2,3,  and  4).  The  over¬ 
shoot  is  little  affected  due  to  A  A^^  and  ^  A^q  and  strongly 
affected  due  to  AwA^o  ^4d  ^  ^Ao* 

The  plot  of  X(,j=  j=1,2,3,  and  4  in  Fig. 3. 7 

indicates  that  X6(t)  is  little  affected  in  the  rise  time  due 
to  parameter  variations  A  A j  ^  (j=1,2,3,  and  4).  The  over¬ 
shoot  and  steady  state  are  little  affected  due  to  A  A 


(O 


and 


A  A 2o  strongly  affected  due  to  A  a^q  ^nd  A  A^^. 

Table  I  shows  the  above  analysis  in  a  condensed  way  to 
give  a  general  picture  of  all  states  and  sensitivity  func¬ 
tions  with  the  correspondent  effect  as  a  function  of  time. 


E.  PARAMETEH-IHDOCED  OOTPOT  ANALYSIS 

If  A«<  «  Xo  ,the  actual  output  can  be  written  as 
(Eqn.  2.29)  . 


y(t,X)  =  y(t,o<^) 


♦  ^ (t,o^)  AX 


(3.17) 


The  computer  program  as  shown  in  the  Appendix  D  gives 
the  actual  output  when  10%,  3015  and  40%  of  variation  from  the 
nominal  value  of  each  parameter  is  assumed.  Figs. 3.8 
through  3.13  show  the  actual  output  when  each  parameter  is 
varied  10%  from  the  nominal  value. 

Fig. 3. 8  indicates  that  the  overshoot  and  rise  time  of 
the  actual  and  nominal  output  of  XI  are  strongly  affected 
which  is  in  agreement  with  the  previous  analysis. 
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?ig.3.9  iadicates  that  the  steady  state  and  overshoot  of 
the  actual  and  nominal  output  of  X2  are  primarily  affected 
which  is  in  agreement  with  the  previous  analysis. 

Fig. 3. 10  indicates  that  the  staady  state  and  overshoot 
of  the  actual  and  nominal  output  of  X3  are  primarily 
affected  which  is  in  agreement  with  the  previous  analysis. 

Fig. 3. 11  indicates  that  the  overshoot  of  tha  actual  and 
nominal  output  of  X4  is  primarily  affected  which  is  in 
agreement  with  the  previous  analysis. 

Fig. 3. 12  indicates  that  overshoot  of  the  actual  and 
nominal  output  of  X5  are  primarily  affected  which  is  in 
agreement  with  the  previous  analysis. 

Fig. 3. 13  indicates  that  the  steady  state  and  overshoot 
of  the  actual  and  nominal  output  of  X6  are  primarily 
affected  , having  little  effect  in  the  rise  time  which  is  in 
agreement  with  the  previous  analysis. 

Figs. 3. 14  through  3.19  and  Figs. 3. 20  through  3.25  show 
respectively  the  actual  output  for  30%  and  40%  of  the 

nominal  value. 

From  the  plots  one  can  see  that  for  small  parameter 
variations  the  parameter-induced  output  error  is  negligible 
and  when  the  variation  becomes  large  as  30%  or  40%  one  note 
that  the  error  becomes  pronounced  and  that  modeling  is 
starting  to  break  down.  This  agrees  with  the  assumption 
made  in  the  derivation  of  the  Eqn.2.28. 

P.  COOPLED  ROLL-iaa  AUTOPILOT  ANALISIS 

1.  Linear  Equations  of  the  Nominal  System 

From  the  block  diagram  of  the  Fig.B.1  in  the 

appendix  B  one  can  get  the  following  state  equations  of  the 

coupled  roll-yaw  autopilot. 
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XI  =  B  conv  (A4  X3  +  A3  X8  ♦  A5  XII) 


(3. 18) 


X2  =  C  Conv  (A7  X3  +  A8  X8  ♦  A6  XII)  (3.19) 

X3  =  -  XI  -  (ALPHAB/Conv)  X2  (3.20) 

+  KB  A  A2  X3  +  KB  A  A1  XII 

X4  =  -  8  X4  -  17.6  X12  +  17.6  PHC  (3.21) 

X5  =  -  (0.755/Coqv)  X2  -  C  D  A7  X3  (3.22) 

♦  (0.755  -  8  D)  X4  -  5  X5  -  C  D  A8  X8  -  C  D  A5  XI  1 

-  17.6  X12  +17.6  PHC 

X6  =  C  E  A7  X3  -  6  X6  ♦  B  C  A8  X8  +  C  E  A6  XII  (3.23) 

« 

X7  =  -  F  KC  (0.755/Conv)  X2  (3.24) 

-  (D  +  E)  F  KC  C  A7  X3  +  P  KC  (0.755  -  8  D)  X4 
+  KC(15  -  5  F)  X5  +  KC(6F  -15)  X6  -  15  X7 

-  (D  +  E)  F  KC  C  A8  X8  -  (D  +  E)  F  KC  C  A6  XII 

-  F  KC  D  17.6  X12  >  F  KC  D  17.6  PHC 

X8  =  188.4  Conv  X7  -  188.4  X8  (3.25) 
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*  »x* 


X9  =  (K1  k  A2/t1)  XI  -  X9  /  T1  +  (K1  k  A1/T1)  XII  (3.26) 


X10  =-  (K2/Conv)  XI  -  (K2  H  ALPHAB/Conv  X2  (3.27) 

(k2/10)  (K1  A  A2/T1)  -»•  B  A4  -  H  ALPHAB  C  A7  X3 
♦  (K2/10)  (B  A3  -  H  ALPHAB  C  A8  X8  +  K2(1 

-  1/(10  T1))  X9  +  (K2/10)  (  (K1  A  A1) /T1  +  B  A5 

-  H  ALPHAB  C  A6)  XII 

XII  =  188.4  Conv  X10  -  188.4  XII  (3.28) 


X12  =  X2  /  Conv  (3.29) 


Eqn.3. 30  gives  the  matrix  representation  of  the 
state  variables  above  mentioned. 


(3.30) 
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The  corr espondence  of  the  states  is; 

X1  =  r,X2  =  p,X3  =  fi  , 

X4  =  X  ,  X5  =  T1  ,  X6  =  XI  , 

X7  =  ^R.c  ,  X8  =  6«.  ,  X9  =  Y  , 

X10  =  Syc  '  X11  =  Sy  ,  X12  =  4> 

The  definition  of  the  constants  Cl  throagh  C25  are 
given  in  appendix  B,  The  parameters  of  interest  for  this 
system  are  given  by  A1  through  a8  which  are  : 


A1  =  Cv,  ,  A2  =  Cv  ,  A3  =  ,  kH  =  C  ,  kS 

A6  =  C-  ,  A7  =  C,„  ,  A8  =  Co  . 

2.  Sensitivity  Equations 


C 

n  6y 


f 


As  showed  in  the  previous  analysis  here  one  can 
apply  the  same  procedure  using  the  TRAJECTORY  SENSITIVITY 
EQUATION  as  given  in  Equ.2.  13. 


“  — —  7)a  —  —  — 

V-:--LXo  --C - (3.31) 

For  the  purpose  of  this  procedure  again  one  can  see 
in  Eqn.3.30  the  correspondent  matri£es  A^,  and  B^.  ^ 

The  partial  derivatives  and  the  are 

evaluated  considering  the  parameter  A  of  interest.  In  this 
case  they  are  A1 ,A2 , . .. , A8,  that  are  respectively  the  aero¬ 
dynamic  coefficients  present  in  the  coupled  roll-yaw  auto¬ 
pilot. 

Once  again, applying  the  partial  derivatives  with 
respect  to  the  parameters  one  can  obtain  eight  matrices 
respectively  as  found  , similarly,  in  the  previous  case  of 
the  uncoupled  pitch  autopilot. 

Here, one  can  see  that  xr — r  =  0  because  B  is  inde- 

o 

pendent  of  the  parameters  of  interest. 

In  terms  of  components  one  obtains 
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(3. 32) 


i=1 ,2, . 12 

j=1,2, . 8 

Here, one  can  see  that  eight  models  are  necessary  to 

stady  the  effect  of  parameter  variations.  These  models  are 

shown  in  Fig. 3. 27.  For  one  parameter  of  interest,  say  A1, 

one  have  the  following  sensitivity  equations  : 

'f»Tr  TP  ir  nr,‘ 

A. 


11 
Azi 

A31 
A41 
A51 
Xfei 

Vi 

Afti 

* 

A^9I 

'^10i 

^111 

Aizf 


Ml 


e 


\fi  "*Ai 


1 

X6 

X7 

^6 

Xg 

X,c 

X11 

Alt 


(3.33) 


’III  '”^1X0 

similar  procedure  can  be  male  for  the  others  parame¬ 
ters  A2  through  A8. 


G.  SENSITIVITY  ANALYSIS 

In  order  to  simulate  simultaneouly  the  nominal  and  the 
sensitivity  models, a  computer  program  was  written  as  shown 
in  the  appendix  E. 

For  analysis  purpose, each  parameter  was  varied  simulta¬ 
neously  10%  from  the  nominal  value. 
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1. 


Analysis  of  o< -Parameter  Variations 

^  X, 


The  plots  of  A,,’  =  —  Ia'  »  j=1,2,.,.,8  in 

-*  o  Ai  o 

Fig.  3.27  ana  3.28  indicate  that  a  parameter  change 
(j=  1 , 2 ,. ..  ,  8)  primarily  affect  the  overshoot  of  1C  1  (t)  . 

The  plots  of  A-i  i'  =  Ia."  »  j=1»2,...,8 


in 


-hAj 

Fig.  3.29  and  3.30  indicate  that  a  parameter  change  AAJf,, 

(j=1 , 2,. .. ,  8)  primarily  affect  the  overshoot  of  X2(t). 

The  plots  of  A^j-  =  j=1,2,...,8  in 

Fig.  3.31  and  3.32  indicate  that  a  parameter  change  , 

(j=  1 , 2,. . . ,  8)  primarily  affect  the  overshoot  of  X3(t). 

The  plots  of  A^j  =  j=1,2,...,8  in 

Fig.  3.33  and  3.34  indicate  that  a  parameter  change 

(j=1  ,2  ,. . . ,  8)  primarily  affect  the  overshoot  of  X4(t). 

The  plots  of  Agj  =  j=1r2,...,8  in 

Fig. 3.35  and  3.36  indicate  that  a  parameter  change  A 

(j=1 , 2 ,. . . ,  8)  primarily  affect  ;!:he  overshoot  of  X5(t). 

X 

The  plots  of  o'  j=1»2,...,8  in 

Fig.  3.37  and  3.38  indicate  that  a  parameter  change 

(j=1 , 2 ,. . . ,  8)  primarily  affect  the  overshoot  of  X6  (t)  . 

The  plots  of  Ayj’  =  j=1r2,...,8  in 

Fig.  3.39  and  3.40  indicate  that  i  parameter  change 

(j=  1 »  2 ,. . . ,  8)  primarily  affect  the  overshoot  of  X7(t). 

The  plots  of  /Afij’  =  j=1,2,...,8  in 

Pig.  3.41  and  3.42  indicate  that  a  parameter  change 

(j  =  1 » 2  ,. . . ,  8)  primarily  affect  the  overshoot  of  X8(t). 

The  plots  of  AgJ  =  j=1,2,...,8  in 

Fig. 3.43  and  3.44  indicate  that  a  parameter  change  A  AJq, 

(j=1 » 2 ,. .. , 8)  primarily  affect  the  overshoot  and  rise  time 

of  X9(t)  .  ^  ^ 

I  U  X  j  o 

The  plots  of  A|oj’  =  j=1,2,...,8  in 

Fig.  3.45  and  3.46  indicate  that  a  parameter  change 
(j=1 »  2,. .. ,  8)  primarily  affect  the  rise  time  of  Xl0(t). 
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The  plots  of  -^1 


j=l,2,...,8  in 


^  •  J  ^  A  j  Aj  o  ' 

Fig.  3.47  and  3.48  indicate  that  a  parameter  change 

(3=“* » 2  ,.  . . , 8)  primarily  affect^he  overshoot  of  lC11(t). 

The  plots  of  /\  ( 2^  —  “  I  /\  j  3“^»2»...r8  in 

Fig.  3.49  and  3.50  indicate  that  a  parameter  change 
(j=1 , 2 ,  8)  primarily  affect  the  overshoot  of  X12(t). 

Table  II  and  III  show  the  above  analysis  in  a 
condensed  way  to  give  a  general  picture  of  all  states  and 
output  sensitivity  functions  with  the  correspondent  effect 
as  function  of  time. 

2.  Parameter- Induced  Output  Analysis 

As  Shown  previously,  if  ^  oc  <<  ,  the  actual 

output  can  be  written  as 


y(t,o<.)  =  y(t, 


+  ^  (t,  c^)  ^ 


(3.34) 


The  computer  program  given  in  appendix  E  was  written 
for  simulating  the  system  when  10??,  30%,  and  43%  of  varia¬ 

tion  from  the  nominal  value  of  each  parameter  is 
assumed. Fig. 3.  51  through  3.62  give  the  actual  output  when 
each  parameter  is  varied  10%  from  the  nominal  value. 

Fig. 3. 63  and  3.64  give  the  actual  output  for  30%  of 
variation  from  the  nominal  value  of  each  parameter  assumed. 
Fig. 3.65  and  3.66  give  the  actual  output  for  43%  of  varia¬ 
tion  from  the  nominal  value  of  each  parameter  assumed. 
These  plots  show  the  output  of  the  state  variables  X3  and 
XII  that  present  strong  deviations  just  to  give  the  behavior 
of  the  system  when  parameter  variations  are  not  small.  One 
notes  that  modeling  is  starting  to  break  down. 
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Figure  3,2  Sensitivity  of  XI  with  Hespect  to  A1,&2^A3,A4. 
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Figure  3.3  Sensitivity  of  12  with  Respect  to  ai,h2,A3,A4. 
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Figure  3,4  Sensitivity  of  X3  with  Respect  to  &1,A2,A3,A4. 
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Figure  3.5  Sensitivity  of  X4  with  Respect  to  A1,A2,A3,A4. 
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Figure  3.6  Sensitivity  of  X5  with  Respect  to  A1,A2,A3,A4. 
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Figure  3.8 


Actual  and  Noainal  Output  of  XI 


(10X  variation)  . 
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Figure  3,9 


ictual  and  Nominal  Output  of  12  (10%  variation) 


TIME 


I 

I 


LEGEND 


Figure  3.10 


Actual  and  Nominal  Output  of  X3 


(10%  variation)  . 
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Figure  3.11 


Actual  and  Noiinal  Output  of  X4 


(10%  variation)  . 
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Figure  3.12 


Actual  and  Nominal  Output  of  X5 


(10?5  variation)  . 
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13  Actual  and  Nominal  Output  of  X6  (10%  variation) 
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Figure  3.14  Actual  and  Nominal  Output  of  XI  (33%  variation). 
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Figure  3.15  Actual  and  Noninal  output  o£  Z2  (33%  variation) . 
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Figure  3.16 


Actual  and  Nominal  Output  of  X3 


(30  X  variation)  . 
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Figure  3.17 


Actual  and  Nominal  Output  of  Z4 


(30 variation)  . 
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Figure  3.18 


Actual  and  Noainal  Output  of  15 


(30 X  variation)  . 
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Figure  3.19 


Actual  and  Nominal  Output  of  Z6 


(3^%  variation)  . 
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Figure  3^20  Actual  and  Noainal  Output  of  XI  (43%  variation) 
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Figure  3.21 


actual  and  Nominal  Output  of  X2 


(49%  variation)  . 
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Figure  3.22 


Actual  and  Nominal  Output  of  X3 


(43  36  variation)  . 
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23  Actual  and  Noninal  Output  of  X4  (40%  variation) . 
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Figure  3.24 


Actual  and  Nominal  Output  o£  X5 


(49 X  variation)  . 
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3,25  Actual  and  Nominal  Output  of  X6  (43  56  variation) 
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Figure  3.26  Boll-Yav  Nominal  and  Sensitivity  Hodels. 
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3.27  Sensitivity  of  XI  with  Respect  to  h1,A2,A3, 44 
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Figure  3.28  Sensitivity  of  XI  with  Respect  to  i5,A6,A7, A8. 
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Figure  3.29  Sensitivity  of  X2  with  Respect  to  &1rA2, 43,14. 
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Figure  3.30  Sensitivity  of  X2  with  Respect  to  &5,46,A7,A8. 
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Figure  3.31  Sensitivity  of  X3  with  Respect  to  &1,A2,A3,A4. 
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Figure  3.32  Sensitivity  of  13  witli  Respect  to  &5,A6,A7,  A8. 
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3.33  Sensitivity  of  X4  with  Respect  to  M,A2,A3,A4. 
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Figure  3.34  Sensitivity  of  X4  with  Respect  to  45,A6,a7,A8. 
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Figure  3.35  Sensitivity  of  X5  with  Respect  to  M,A2,A3,A4. 
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Figure  3.36  Sensitiyity  of  X5  with  Respect  to  h5,A6,A7, AS. 
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Figure  3,37  Sensitivity  of  X6  with  Respect  to  hi ,12,43, A4. 
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Figure  3.38  Sensitivity  of  16  with  Respect  to  i5,A6,A7,A8. 
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Figure  3.39  Sensitivity  of  X7  with  Bespect  to  &1«A2,A3,A4. 
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Figure  3.40  Sensitivity  of  17  witli  Respect  to  A6,  A7,  AS. 
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Figure  3,41  sensitivity  of  X8  with  Respect  to  ai,a2,A3eA4. 
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Figure  3*42  Sensitivity  of  X8  witli  Respect  to  i5,A6,A7,A8, 
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Figure  3.43  Sensitivity  of  X9  vitti  Respect  to  &1,A2rA3,A4. 
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Figure  3.44  Sensitivity  of  X9  with  Respect  to  &5,A6,A7,A8. 
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Figure  3.45  Sensitivity  of  X10  with  Respect  to  i1,A2,A3,a4. 
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Figure  3.46  Sensitivity  of  X10  with  Bespect  to  A5,a6,a7,A8. 
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Figure  3.47  Sensitivity  of  XII  with  Respect  to  A1,A2,A3,A4 
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3.48  Sensitivity  of  X11  with  Respect  to  A5,A6,A7,A8, 
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Figure  3.49  Sensitivity  of  X12  with  Respect  to  A1rA2rA3^A4, 
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Pigure  3.50  Sensitivity  of  X12  with  Respect  to  a5,A6,A7,A8. 
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3*51  Actual  and  Nominal  Output  of  XI  (10J6  variation) 
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Figure  3.52  Actual  and  Nominal  Output  of  12^^0%  variation). 
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Pigure  3.53  Actual  and  Nominal  Output  of  X3(10*  variation). 
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Figure  3.54  Actual  and  Nominal  Output  of  X4  (10X  variation)  . 
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Figure  3*55  Actual  and  Nominal  Output  of  X5(10X  variation). 


109 


4 


LEGEND 


110 


i 


I 

i 


60*0  iO‘0  90 '0  fiO'O  £0‘0  TpO  00*0-  20*0-  £0*0- 


Figure  3.57  Actual  and  Nominal  Output  of  X7(10S  variation). 
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Figure  3*58  Actual  and  Nominal  Output  of  Xd(10%  variation). 
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Figure  3*59  Actual  and  Nominal  Output  of  X9  (10X  variation). 
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Figure  3.60  actual  and  Nominal  Output  of  X10(13*  variation). 
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Figure  3.62  Actual  and  Hominal  Output  of  X12(1(}%  Trariation)  . 
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Figure  3.63  Actual  and  Nominal  Output  of  13  (30it  Tariation)  . 
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Figure  3.64  Actual  and  Noninal  Output  of  XI 1(30%  variation). 
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Figure  3*65  Actual  and  Nominal  Output  of  X3  (40S  variation). 
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Figure  3.66  Actual  and  Nominal  Output  of  111(43%  variation). 
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17.  ftPPLICATIOH  OF  SENSITIVITY  ANAtXSIS  TO  NO^INEAR  SYSTEMS 

A.  INTRODDCTION 

The  three  dimensional  nonlinear  system  used  here  for 
sensitivity  analysis  purposes  is  presented  in  appendix  c. 
The  missile  configuration,  size  and  mass  properties  are  the 
same  as  those  used  in  the  lineat  system  and  were  presented 
in  appendix  A, 

In  section  B  the  sensitivity  equations  with  respect  to 
parameters  of  nominal  system  is  shown  and  in  section  C  the 
sensitivity  equations  with  respect  to  parameters  that  showed 
to  be  of  more  effect  in  the  previous  analysis  is  presented. 

Section  D  shows  the  results  of  the  trajectory  sensi¬ 
tivity  equations  when  step  inputs  are  applied  at  specific 
time  as  mentioned  in  appendix  C. 

Section  E  gives  the  parameter-induced  output  analysis 
when  each  parameter  is  varied  10%  from  the  nominal  value. 

B.  NONLINEAR  EQUATIONS  OF  THE  NOMINAL  SYSTEM 

From  the  block  diagram  of  Fig.  C.1  and  Egn.  C.3  through 

C. 21  as  given  in  appendix  C  one  have  the  following  nominal 
equations. 


XI  =  -  C4  XI  -  C02  C4  A2  (4.1) 


X2  =  -  C5  X2  -  C6  NZC  ♦  C6  C19  Cos  (X7)  +  C7  XI  (4.2) 
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X3  =  (C5  C8  -  C9)  X2  +  C6  C8  NZC  (4.3) 

-  C6  C8  C19  Cos(X7)  -  C7  C8  XI  +  (C8/C0liv2)  XI  7  X18 
+  Cl  C8  A1  +  (C9/Conv)  X5 

X4  =  -  C3  X4  +  C3  Conv  X3  (4.4) 

X5  =  (X17  X18)/Conv  +  Cl  Conv  A1  (4.5) 

X6  =  X5  -  KB  C02  A2  -  (X16  X18)/Conv  (4.6) 

X7  =  X5  Cos (XI 9) /Conv  -  X17  Sin  (XI  9) /Conv  (4.7) 

# 

X8  =  -  CIO  X8  ♦  (C12  -  Cl  1  C13)  XIO  (4.8) 

+  C11  C14  PHC  -  C11  C14  X19  -  C11  CO (A8  X16  +  A7  X15 
♦  A3  X13)  -  (C12/Conv)  X18 

X9  =  -  X9/T1  +  K1/T1  C02(AA  X16  ♦  A6  X15)  (4.9) 

4 

XIO  =  -  C13  XIO  +  C14  PHC  -  C14  X19  (4.10) 

XII  =  -  C15  XII  +  C16  C0(A8  X16  ♦  A7  X15  +  A3  X  13)  (4.11) 


122 


X12  =  -  C17  X12  +  K(C17  -  (C10/C18)  X8 


(4.  12) 


(K/C18)  (C12  -  C11  C13)  X10  ♦  K  Cl  1  (C14/C18)  PHC 

-  K  C11  (C14/C18)  X19  -  (K/C18)  CO  A8(C11  -  (K/C18) 

CO  A3(C11  *■  C16)  X13  -  (K/C12)  (Cl  8  Conv)  X18 

+  K(C15/C18)  -  C17)  XII 

X13  =  -  C3  X13  ♦  C3  Conv  X12  (4.13) 


X14  =  (K2/10)  (-  (X9/T1)  C02(AA  X16  ♦  A6  X15)  (4.14) 

-  X15  (X18/Conv2)  +  C01(A9  X16  ♦  A4  X15  ♦  A5  X13) 

-  (KYP/Conv2)  (X5  -  KB  C02  A2  -  XI 6  (X 1 8/Con v)  X18 

-  (KYP/Conv)  XX6  C0(a8  X16  ♦  A7  X15  +  A3  X13) 

♦  (K2/Conv)  X17  -  K2(KYP/Conv2)  XX6  X18 

X15  =  -  C3  X15  ♦  C3  Conv  X14  (^.15) 


X16  =  KB  C02(AA  X16  ♦  A6  X15)  ♦  X5(X18/Conv)  -  X17(4.16) 

X17  =  -  (X5  X18)/Conv  +  C01  Conv(A9  X16  (4.17) 

+  A4  X15  +  A5  X13) 

X18  =  CO  Coav(A8  X16  +  A7  X15  A3  X13)  (4.18) 


XI 9  =  Xl8/Conv 


(4.  19) 


The  correspondence  of  nonlinear  state  vectocsi  as  given 
in  the  above  equations  with  relation  to  the  nonlinear  system 
presented  in  Fig.C.  1  in  appendix  c  are  : 

XI  =  X  ,  X2  =  Y  ,  X3  =  ,  X4  =  Sp  ,  X5  =  q  , 

X6  =  ,  X7  =  ©  ,  X8  =  Y1  ,  X9  =  Y2  ,  XI  0  =  XI  , 

XII  =  X2  ,  X12  =  Src»  X13  =  Sr,  X14  =  X15  =  8^, 

X16  =  ,  X17  =  r  ,  X18  =  p  ,  and  X19  =  4  . 


The  parameters  of  interest  for  this  nonlinear  system  are 
given  by 

A1  a2  =C}<X,ip\,  A3  =C„U),  a4  =  ,  and  A5 


M 


"«Y 

Definition  of  the  constants  CO  through  C19  are  given  in 
appendix  C. 

As  previously,  this  nomenclature  is  used  here  to  easily 
apply  the  sensitivity  theory. 


C.  HOHLIHBAR  SENSITIVITY  EQUATIONS 

From  the  sensitivity  theory  given  in  chapter  2  one  knows 
from  Egn.2.12  that  for  nonlinear  systems  one  has 

'1)X 

^ - =  — - + - -  r - =  0  (4.20) 

Here 


X  =  f  {X1,X2,..  .,X19,t»NZC,Al,A2,...  ,A5} 


(4.  21) 


iThe  state  vectors  and  parameters  in  the  nonlinear  case 
have  no  correspondence  to  the  stats  vectors  and  parameters 
used  in  the  linear  case. The  correspondence  can  be  found  in 
the  appendices  B  and  C. 
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From  Eqn.2.14  one  knows  that  the  trajectory  sensitivity 
vector  is  given  as 


.li 


C<0 


A  (0)  =  0 


(4.22) 


Applying  the  above  theory  one  obtains  the  nonlinear 
sensitivity  equations  as  given  in  Eqn.4.24  through  54942, 
where  ^  was  considered  the  parameter  to  be  varied. 


Ill  =  -  C4/111  -  C02  C4{DA26a61  ♦  DA24a41) 


(4. 23) 


A21  =  -  C5/i21  -  C6  C19  Sin(X7)  A71  +  C7^11  (4.24) 


,131  =  (C5  C8  -  C9)  ,121  +  C6  C8  C19  Sin  (X7)  ,171  (4.25) 

-  C7  C8  All  +  C8(x18A171  (X17  A  18  1 )  /  (ConvZ) 

+  Cl  C8(DA16A61  +0114  A41  +  1.)  +  (C9A51)/Conv 

A41  =  -  C3A41  +  C3  Conv  A31  (4.26) 


A51  =  (X17A181  ♦  X18^17 1)/Conv  (4.27) 

+  Cl  Conv(DA16A61  +  DA14A41  +  1.) 

A61  =A51  -  KB  C02(DA26/i61  +  DA24A41)  (4.28) 

-  ((X16A181  +  X18  A161) /Conv 
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.171  =  (Cos(X19)  /i51  -  X5  Sin  (X19)  .1191) /COQV  (4.29) 

-  (Sin(Xl9)  .1171  +  X17  Cos  (X19)  ai71) /conv 

Jlai  =  -  CIO  1.81  +  (C12  -  C11  C13)  A101  (4.30) 

-  C11  C141.191  -  C11  CO  (A8/i161  ♦  A7.i151  +  43^131 

+  /d61  (X16  Da8  X15  Da7  4X13  DA13))  -  (Cl  2  1.1 81) /Conv 

1.91  =  -1L91/T1  ♦  K1/(T1  C02)  (AAl.161  (4.31) 

+  A6  1.151  *■  1.61  (X16  DAA  +  X15  Da6)  ) 

# 

1.101  =  -  C13  1.101  -  C141.191  (4.32) 

/ 

m 

1111  =  -  C15  1.111  +  C16  C0(A81161  +  A7H51  (4.33) 

+  A31131  ♦  161  (X16  DA8  ♦  X15  DA7  +  X13  DA3)  ) 

1121  =  -  C17  1121  +  K(C17  -  C10/C18)  181  (4.34) 

+  (K/C18)  (C12  -  C11  C13)  1101  -  K  C11  (C14/C18)  1191 

-  (K/C18)  C0(C11  ♦  C16)  (A81161  +  A7  1151  +  A3  1131 
♦  A3  1131  +  ^61(X16  DA8  ♦  X15  DA7  +  X13  Da3)  ) 

-  K  C12/(C18  Conv)  1181  +  K  (C15/C18  -  C17)1111 

• 

1131  =  -  C31131  +  C3  ConvH21  (4.35) 
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^141  =  (K2/10)  (-  X9/T1  +  (K1/T1)  C02(AA 


(4.36) 


/i161  +  A6/il51  +  261(X16  DAA  +  X15  DA6)  )  -  (X15/i181 
+  X18  Z151)  /  (C0I1V2)  +  C01(A9.1161  +  A4^151  +  AS 
+  ^61(X16DA9  +  X15  DA4  +  X13  DA5)  )  -  KYP/(Coav2) 

(X5i.181  ♦  X18/i51  -  KB  C02(A2^181  +  X18(DA25^61 
♦  DA24/i41)  -{X182ai61  +  2  X18  X16/1181)  /Conv)) 

-  (KYP/Con\r)  C0{  xx6/i161  +  X16ai61)  +  A7ai61) 

+  A7(XX6/1151  ♦  Xl5/il61)  +  A3(XX6/i131  +  X13/1161) 

+  X16(XX6  DA8/i61  +  A8^161)  +  X15{XX6  DA7xl61 
+  A7  Z161)  +  X13(XX6  DA3  ^61  +  A3  ^161))) 

+  (K2/Conv)  /1 171  +  K2/i91  -  K2  KYP/(Conv2) 

(XX6  /1181  ♦  X18  ^161)  ) 

^151  =  -  C3yl151  +  C3  Coavi.141  (4.37) 


A161  =  KB  C02(AA  ^61  +  A6  2.151  +^61(  (4.38) 

X16  DAA  +  X15  DA6))  +  (X6/1181  +  X 1 8  /16 1) /Conv  -^171 

/1171  =  -  (X5  2.181  +  X18  2.51) /Conv  (4.39) 

+  C01  Conv(A92.161  +  A4  2.151  +  A5xi131  +2.61(X16  DA9 
+  X15  DA4  ♦  X13  DA5)  ) 


2.181  =  CO  conv  (A8  2.161  +  A7  2151 


(4.40) 


+  A32.131  +2L61(X16  DA8  ♦  X15  DA7  +  X13  DA3)) 


2191  =2181/Conv 


(4. 41) 
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The  same  procedure  above  must  be  made  for  the  other 
parameters  that  most  affected  the  previous  analysis  and  will 
be  given  in  the  nonlinear  sensitivity  analysis  presented  in 
the  next  section. 


the  nominal 
Here,  the 
cases  given 
in  the  time 


D.  NONLIHEAR  SENSITIVITY  ANALYSIS 

A  computer  program  in  appendix  3  shows  the  simulation  of 
the  sensitivity  eguations  with  respect  to  the  selected 
parameters. 

The  number  of  equations  solved  are  19  for 
equations  and  95  for  the  sensitivity  equations, 
parameters  were  selected  from  the  previous  two 
in  Chapter  III  that  showed  to  have  most  effect 
response  of  the  systems. 

From  the  analysis  performed,  two  paraueters  were 
selected  from  the  uncoupled  pitch  autopilot  and  three  from 
the  coupled  roll-yaw  autopilot,  which  has  been  shown  to  be 
more  sensible  to  parameter  variations.  These  parameters  are 
respectively  (A2)  and  (A4)  in  the  uncoupled  pitch 

autopilot  and  ( A3)  ,  (^5)  ,  and  in  the 

coupled  roll-yaw  autopilot. 

In  order  to  compare  the  influence  of  the  parameter  of 
interest  in  the  nonlinear  case  some  state  variables  were 
chosen  to  be  analysed.  The  state  variables  selected  were: 

Prom  the  uncoupled  pitch  autopilot  : 

X3  =  Spg  ,  X4  =  Sp  ,  X5  =  q  ,  X6  =  X 

From  the  coupled  roll-yaw  autopilot  : 

X12  =  <Sr^  ,  X13  =  ,  X14  =  ,  X15  =  Sy  , 

X16  =  ,  X17  =  r  ,  X18  =  p  ,  and  X19  =  . 
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The  results  are  plotted  in  Fig. 4.1  through t  4.12  and 
Table  IV  and  V.  Each  state  variable  output  is  plotted 
separately  with  the  correspondent  five  output  sensitivity 
functions. 

From  the  plots,  the  following  observations  can  be  made: 

Fig. 4.1  indicates  that  the  state  variable  X3(5f>^)  is 

little  affected  in  the  rise  time  by  the  parameters  A1  and 
A2.  The  overshoot  is  strongly  affected  by  A5,  A1  and  A4  and 
little  affected  by  A2  and  A3.  The  steady  state  is  strongly 
affected  by  A1  and  little  affected  by  A2. 

Fig.  4. 2  indicates  that  the  state  variable  X4(<Sp)  is 

little  affected  in  the  rise  time  by  A1  and  A2.  The  over¬ 
shoot  is  strongly  affected  by  A1,  A3,  and  A4  with  little 

effect  by  A2  and  A5 .  The  steady  state  is  strongly  affected 
by  A1  and  little  affected  by  A2. 

Fig. 4. 3  indicates  that  the  state  variable  X5  (  cj  )  is 

strongly  affected  in  the  rise  time  by  A1  and  A2.  The  over¬ 
shoot  is  strongly  affected  by  A3,  A4,  A5,  and  little 

affected  by  A1  and  A2.  The  steady  state  is  not  affected  by 

parameter  variations. 

Fig. 4. 4  indicates  that  the  state  variable  X6  (  )  is 

strongly  affected  in  the  rise  time  by  A1  and  A2.  The  over¬ 
shoot  is  strongly  affected  by  A3,  A4,  and  A5.  The  steady 

state  is  little  affected  by  A1  and  A2. 

Fig. 4. 5  indicates  that  the  rise  time  and  steady  state  of 
the  state  variable  X12(5p^J  are  not  affected  by  parameter 
variations.  The  overshoot  is  strongly  affected  by  A3  and  A5 
and  little  affected  by  A1,  A2,  and  A4. 

Fig. 4.6  indicates  that  the  rise  time  and  steady  state  of 
the  state  variable  X13(Jrj^)  are  not  affected  by  parameter 
variations.  The  overshoot  is  strongly  affected  by  A3,  A4, 

and  A5  and  little  affected  by  A1  and  A2. 

Fig.  4. 7  indicates  that  the  rise  time  and  the  steady 
state  of  the  state  variable  X14(Sy^)  are  not  affected  by  the 
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parameter  variations.  The  overshoot  is  strongly  affected  by 
A4  and  A5  and  little  affected  by  A1,  A2,  and  A3. 

Fig. 4. 8  indicates  that  the  rise  time  and  the  steady 

state  of  the  state  variable  X15(6y  >  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  A5  and  little  affected  by  Al,  and  A2. 

Fig. 4. 9  indicates  that  the  rise  time  and  the  steady 

state  of  the  state  variable  X16(^  )  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  A5  and  little  affected  by  Al,  and  A2. 

Fig. 4. 10  indicates  that  the  rise  time  and  the  steady 

state  of  the  state  variable  X17(  )  are  not  affected  by  the 

parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  A5  and  little  affected  by  Al,  and  A2. 

Fig. 4. 11  indicates  that  the  rise  time  and  the  steady 

state  of  the  state  variable  X18(  \  are  not  affected  by  the 

parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  A5  and  little  affected  by  Al,  and  A2. 

Fig. 4. 12  indicates  that  the  rise  time  and  the  steady 

state  of  the  state  variable  X19(  (j>  |  are  not  affected  by  the 
parameter  variations.  The  overshoot  is  strongly  affected  by 
A3,  A4,  and  A5  and  little  affected  by  Al,  and  A2. 

E.  PABAHETER-INDOCED  OOTPOT  ANALYSIS 

As  shown  previousy,  if  o^q  ,  the  actual  output  can 

be  written  as 


y  t  y  (t,  )  *  (T  {t ,  o<^  )  ac 


(4.42) 


The  computer  program  given  in  appendix  G  was  written  for 
simulating  the  system  when  variations  from  the  nominal  value 
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of  each  parameter  is  assumed.  Figs. 4. 13  through t  4.24  give 
the  plots  of  the  actual  and  nominal  output  of  the  state 
variables  that  are  present  in  the  unclouped  pitch  autopilot 
and  coupled  roll-yaw  autopilot  when  10S  of  parameter  varia¬ 
tion  is  assumed.  The  others  state  variables  are  not  given 
here  because  they  are  just  intermediate  state  variables. 

these  results  confirm  the  theory  presented  in  chapter  2 
when  small  parameter  variations  are  assumed.  Figs. 4.25  and 
4.26  present  the  actual  and  nominal  output  of  the  state 
variables  X3  and  XII  that  showed  pronounced  variations  when 
30K  of  parameter  variation  is  assumed.  Figs. 4.27  and  4.28 
present  the  actual  and  nominal  output  of  X3  and  XII  when  40% 
of  parameter  variation  is  assumed.  From  these  polts  one 
note  that  the  modeling  is  starting  to  break  down. 
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Figure  4.1  Sensitivity  of  X3  with  respect  to  &1  , A2«A3,A4, AS. 
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Figure  4.3  Sensitivity  of  X5  with  respect  to  &1 ,&2,A3,A4,A5. 
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Figure  4.4  Sensitivity  of  X6  with  respect  to  41 ,A2,A3,A4, A5- 
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Figure  4.5  sensitivity  of  XI 2  with  respect  to  i 1 , A2, &3, A4, A5. 
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Figure  4.6  Sensitivity  of  X13  with  respect  to  M,A2, A3, A4,i5. 
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Figure  4.7  Sensitivity  of  X14  with  respect  to  M,A2, A3,A4, &5. 
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Figure  4*8  Sensitivity  af  X15  with  respect  to  M ^ A2^ A3, A4^ AS. 
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Figure  4,9  Sensitivity  of  X16  with  respect  to  &  1 A3, A4, &5. 
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Figure  4.10  Sensitivity  of  X17  with  respect  to  A1,A2r A3,A4, AS. 
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Sensitivity  of  X18  with  respect  to  Al, A2, A3,A4^ A5. 
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Sensitivity  of  X19  with  respect  to  Al» A2,A3»i4, A5. 
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Fi.9ur€  Actusil  and  Noni.na.1  Output  of  X3  (13%  variation)  • 
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Figure  4.14  Actual  and  Nominal  Output  of  X4  (13  55  variation)  . 
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Figure  4.15  actual  and  Hoainal  Output  of  X5  (10 S  variation) . 
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4.16  Actual  and  Noainal  Output  of  X6  (13%  variation). 
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Figure  4.17  Actual  and  Noninal  Output  of  XI 2  (10X  variation). 
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18  Actual  and  Nominal  Output  of  XI 3  (10X  variation)* 
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Figure  4,19  Actual  and  Noainal  Output  of  K14  (1  055  variation), 
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Figure  4.20  Actual  and  Nominal  Output  of  X15  (10%  variation). 
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Figure  4.21  Actual  and  Honinal  Output  of  X16  (10*  variation). 
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Figure  4,22  Actual  and  Nominal  Output  of  117  (10X  variation). 
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Figure  4 


.23  Actual  and  Noainal  Output  of  X18  (10%  variation). 
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4.24  Actual  and  Nominal  Output  of  XI 9  (10%  variation). 
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Figure  4.25  Actual  and  Noninal  Output  of  X6  (30 X  variation)  . 
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»,26  actual  and  Nominal  Output  of  X16  (30?t  variation) 
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Actual  and  Nominal  Output  of  X6  (43%  variation). 
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Figure  4.28  Actual  and  Honinal  Output  of  116  (40%  variation). 
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V,  CONCLOSIONS  AND  RBC3C1MENDATIOHS 


A.  CONCLUSIONS 

The  results  of  the  nonlinear  3-D  parameter  sensitivity 
analysis  presented  in  chapter  IV  verify  rhe  linear  analysis 
of  chapter  III.  These  conclusioas  can  be  quickly  made, 
using  Tables  I,  II,  III  and  IV  that  give  a  brief  review  of 
the  time  response  of  the  linear  and  nonlinear  system. 
Another  means  of  comparison  of  the  linear  and  nonlinear 
analysis  is  to  use  the  figures  that  give  a  precise  visuali¬ 
zation  of  the  time  response. 

The  comparison  of  the  linear  and  nonlinear  analysis  can 
be  done  by  means  of  the  plots  as  follows: 


B.  EFFECT  ON  STATE  VABIABLES  DDE  TO  PARAMETER  VARIATIONS 

1 .  Effect  on  State  Variable  Due  to  Variation  in 

tn 

Figs. 3. 7  and  4.1  (  Ac4  and  )  show  that  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  steady  state. 

Figs.  3. 4  and  4.2  (  and  X4  j  )  show  that  ^f>  is 

little  affected  in  the  rise  time  and  strongly  affected  in 
the  overshoot  and  steady  state. 

Figs. 3. 2  and  4.3  (X,^  andAs”))  show  that  q  is 

strongly  affected  in  the  rise  time,  little  affected  in  the 
overshoot  and  not  affected  in  the  steady  stats. 

Figs. 3. 3  and  4.4  (  A^l  ^  show  that  «<  is 

strongly  affected  in  the  rise  time,  little  affected  in  the 
overshoot  and  steady  state. 
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2. 


Effect  on  State  Variable  Due  to  Variation  in 

C_(ot,iE)  : 

M 


Figs.  3. 7 

and 

(  XfcZ  A32)  show 

t  hat 

^Pc  is 

little 

affected 

in 

the  rise  time,  overshoot 

and 

steady 

stat  e. 

Figs.  3 

.4 

and  4.2  (^92.  )  show 

that 

Sp  is 

little 

affected 

in 

the  rise  time,  overshoot 

and 

steady 

state. 

Figs. 3. 2 

and  4.3  (A12  and  show 

that  g  is 

strongly  affected  in  the  rise  time  and  little  affected  in 
the  overshoot  and  steady  state. 

Figs. 3. 3  and  4.4  (^22  >^€*2)  show  that  o<  is 

strongly  affected  in  the  rise  time  and  little  affected  in 
the  overshoot  and  steady  state. 

3.  Effect  on  State  Variable  Due  to  Variatioi  in  C, 

Figs. 3. 40  and  4.6  (A7g  and  Ai23)  show  that  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 42  and  4.7  (Ag^and  A155)  show  that  is  not 
affected  in  the  rise  time  and  steady  state,  and  strongly 
affected  in  the  overshoot. 

Figs. 3. 46  and  4.8  (  Aio©  )  show  that  6yc,  is 

little  affected  in  the  rise  time,  little  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 48  and  4.9  (  Ajig  and  A153  )  show  that  is 

little  effected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs.3.32  and  4.10  (A35andA|63  )  show  that  is 

little  affected  in  the  rise  time,  stronglu  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 28  and  4.11  (A|Q  and  Ai73)  show  that  )r  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 
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Figs. 3. 30  and  4.12  ( Al63  )  show  that  p  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs.  3 . 50  and  4.13  (A(2gaiidAi<^3)  show  that^is 
not  affected  in  the  rise  time  and  steady  stats  and  strongly 
affected  in  the  overshoot. 

4.  Effect  on  State  Variable  Due  ^  Variation  in  C  _  ; 

Figs. 3. 40  and  4.6  and  A|2<q)  show  that  is 

not  affected  in  the  rise  time  and  steady  state  and  little 

affected  in  the  overshoot. 

Figs. 3. 42  and  4.7  (AeS  and  A134  )  show  that  is 

not  affected  in  the  rise  time  and  steady  state  and  strongly 

affected  in  the  overshoot. 

Figs. 3. 46  and  4.8  (XiqS  )  show  that  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  stats. 

Figs. 3. 48  and  4.9  (Aji^  and  A^5^)  show  that  Sy  is 
not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the  overshoot. 

Figs. 3. 32  and  4.10  (  As^and  A  14”^  )  show  that  p)  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  stats. 

Figs. 3. 28  and  4.11  {  A»5  and  Ajy^)  show  taat  r  is  not 

affected  in  the  rise  time  and  steady  srate  and  strongly 
affected  in  the  overshoot. 

Figs. 3. 30  and  4.12  (AzS  and  Aie<^  )  show  that  p  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 50  and  4.13  (A|2,5  Aig^)  show  that  ^  is 

not  affected  in  the  rise  time  and  steady  state  and  strongly 
affected  in  the  overshoot. 
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5.  Effect  on  State  Var iabla  Dus  ^  Variatioi  in  C _ : 

nsp. 

Figs. 3. 39  and  4.6  (^kt^^and  A|2.6‘)  show  that  S(^<;;:,is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 41  and  4.7  (^p^and  A135  )  show  that  is 

little  affected  in  the  risr  time,  stronglu  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 45  and  4.8  (A^o^  ^  show  that  6yc  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 47  and  4.9  (Aji3  and  A/^g’)  show  that  Sy  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 31  and  4.10  (A33ani  A|65  )  show  that  is 

little  affected  in  the  rise  time,  stongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 27  and  4.11  (A13  and  X\qS  )  show  that  V  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 

overshoot  and  not  affected  in  the  steady  state. 

Figs. 3. 29  and  4.12  Ajqs)  show  that  p  is 

little  affected  in  the  rise  time,  strongly  affected  in  the 

overshoot  and  not  affected  in  the  steady  state. 

Figs.3.49  and  4. 13  ( X  123  ^ad  A 13)? )  show  that  ^  is 
little  affected  in  the  rise  time,  strongly  affected  in  the 
overshoot  and  not  affected  in  the  steady  state. 

C.  EFFECT  OF  TRAJECTORY  SENSITIVITY  FONCTIONS 
1 .  Dncoupled  Pitch  Aut opilot 

The  trajectory  sensitivity  functions  ^11,  /i12, 

A13,  and  in  the  linear  case  and  the  correspondent 

trajectory  sensitivity  functions  151  and  152  ia  the  nonli¬ 
near  case  show  strong  effect  on  the  rise  time  and  overshoot 
of  the  state  variable  g. 
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The  trajectory  sensitivity  functions  yi2 3  and  ^24 
in  the  linear  case  and  the  correspondent  trajectory  sensi¬ 
tivity  function  161  in  the  nonlinear  case  show  strong 

effect  on  the  rise  time  of  the  state  variable  pc  . 

The  trajectory  sensitivity  functions  /i33  and  ^34 
in  the  linear  case  and  the  correspondent  trajectory  sensi¬ 
tivity  function  /i41  in  the  nonlinear  case  show  strong 

effect  on  the  overshoot  and  steady  state  of  the  srate  vari¬ 
able  . 

The  trajectory  sensitivity  functions  ^53  and  /i64 
in  the  linear  case  and  the  correspondent  sensitivity  trajec¬ 
tory  function  in  the  nonlinear  case  show  strong  effect 

on  the  overshoot  and  steady  state  of  the  state  variable  • 

2.  Coupled  Ho  11 -Yaw  Autopilot 

The  trajectory  sensitivity  functions  i.13,/i15,  and 
/118  in  the  linear  case  and  the  correspondent  trajectory 
sensitivity  functions /1175,  /11 74,  and/'1173  in  the  nonlinear 

case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  r. 

The  trajectory  sensitivity  functions  ^23 ,  /i25,  and 
2.28  in  the  linear  case  and  the  correspondent  trajectory 

sensitivity  functions  2- 185,  2.184,  and  2.183  in  the  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  p. 

The  trajectory  sensitivity  functions  2.33,  2.35,  and 
238  in  the  linear  case  and  the  correspondent  trajectory 

sensitivity  functions  2165,  >1164,  and2.163  in  the  nonlinear 

case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  ^  . 

The  trajectory  sensitivity  functions  2.73,  275,  and 
278  in  the  linear  case  and  the  correspondent  trajectory 

sensitivity  functions  2 1 25,  2.124,  and  2123  in  tae  nonlinear 
case  show  little  effect  on  the  overshoot  of  the  state  vari¬ 
able 
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The  trajectory  sensitivity  functions  /i83  ,  /185,  and 
^88  in  the  linear  case  and  the  correspondent  trajecrtory 
sensitivity  functions  d.1 35,  /1134,  andxl133  in  tae  nonlinear 
case  show  strong  effect  on  the  overshoot  of  the  state  vari¬ 
able  6(2.  . 

The  trajectory  sensitivity  functions  /II  03  and/1105 
in  the  linear  case  and  the  correspondent  trajectory  sensi¬ 
tivity  functions  /i145  and /II 4 4  in  the  nonlinear  case  show 
strong  effect  on  the  overshoot  of  the  state  variable  • 
The  trajectory  sensitivity  function  /i108  in  the  linear  case 
and  the  correspondent  trajectory  sensitivity  fj net  ion /1 143 
in  the  nonlinear  case  show  little  effect  on  the  overshoot  of 
the  state  variable  Syc.  • 

The  trajectory  sensitivity  functions  /1 113,  /i115, 

and  /1118  in  the  linear  case  and  the  correspondent  trajec¬ 
tory  sensitivity  functions  /1 155,  /1154,  and  /1153  in  the 

nonlinear  case  show  strong  effect  on  the  overshoot  of  the 
state  variable  6y  . 

The  trajectory  sensitivity  functions  /i123,  A125, 

and  /1128  in  the  linear  case  and  the  correspondent  trajec¬ 
tory  sensitivity  functions  /1 195,  /1  194,  and  A193  in  the 

nonlinear  case  show  strong  effect  on  the  overshoot  of  the 
state  variable  . 

The  trajectory  sensitivity  functions  have  the 
following  range  of  values  : 

Minimum  Maximum 

Linear  -  385.25  632.35 

Nonlinear  -  642.9  1772.9 

D.  GENERAL  CONCLUSIONS 

The  parameter  C^^  strongly  affect  the  overshoot  in 
almost  all  case. 
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E.  RECOHBENDiTIOHS 

It  was  not  possible  to  run  a  CSMP  computer  program  using 
all  stare  variable  (19)  and  parameters  (10)  in  the  nonlinear 
case  due  to  the  work  area  available.  The  above  restriction 
occurred  just  when  the  actual  and  nominal  output  were 
required  to  be  printed.  The  computational  time  in  this  case 
may  be  reduced  using  Fortran  subroutines  imbedded  in  the 
CSMP  program  in  order  to  calculate  the  necessary  parameter 
derivatives . 

In  the  Linear  cases  the  computational  time  can  be 
reduced  by  means  of  using  the  method  of  Sensitivity  Points 
that  uses  one  model  for  all  parameters  instead  of  using  one 
model  for  each  parameter. 

The  present  analysis  can  be  repeated  for  the  case  of  the 
Circularz  airframe  given  in  the  [Bef.  2],  and  comparisons 
between  both  airframes  can  be  performed. 

Future  study  can  be  made  modeling  the  system  in  the 
frequency  domain  where  the  "root  sensitivity"  is  analysed. 

Influence  of  parameter  variations  on  miss  distance  can 
be  analysed  using  an  augumented  system  and  a  scenario  of 
reference. 


pr 


2The  current  CSIJP  programs  in  th^  appendices  are 
ared  to  run  the  Circular  case  just  by  inserting  the  data 
he  correspondent  airframe. 
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APPENDIX  A 

MISSILE  SIZING,  HASS  PROPERTIES  AND  AERODYNAHIC  DATA 

A.  INTRODUCTION 

In  this  appendix,  one  considers  a  model  that  was 
assumed  to  be  1/6  scale  of  the  actual  elliptical  missile 
configuration  as  given  in  Fig.A. 1. 

Missile  configuration  was  sized  to  provide  realistc  mass 
properties  needed  for  this  study. 

Since  the  main  purpose  of  this  work  was  to  perform  the 
parameters  sensitivity  analysis  applied  to  a  Sank-ro-Turn 
missile,  no  effort  was  expended  on  a  detailed  design  and 
analysis  of -the  various  parts  envolved,  which  is  given  in 
[Ref.  2].  Some  of  figures  given  in  this  appendix  were  taken 
from  the  [Ref.  2]  fcr  easy  visualization  of  the  system  in 
study. 

B.  GEOMETRIC  AND  MASS  PROPERTIES  OP  MISSILE  CONFIGURATION 

One  can  find  in  Table  VI  the  properties  used  in  the 
development  of  the  equations  for  applying  the  Parameters 
Sensitivity  Analysis. 

Fig. A. 2  shows  the  aerodynamic  sign  convention,  nomencla¬ 
ture  and  body-fixed  axes  where  the  following  three  assump¬ 
tions  were  made: 

-  Plane  x^ -  zl  of  Fig.A. 2  is  the  maneuver  plane. 

-  Missile  is  trimmed  in  pitch  (i.  e.  ,  My=0  ,  at  fixed 

values  of  q,c<,  and  )  . 

-Missile  rcll  rate  is  constant. 

Fig.  A.  3  shows  a  block  diagram  of  a  BTT  autopilot. 
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Figure  a.1  Hodel  of  elliptical  coaf iguration  (  1/6  scale  ). 
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TABLE  VI 

Geoaetric  and  Hass  Properties  of  Missile 
Conf iguratioa. 


WEIGHT  (lbs) .  24  75 

^<3  .  ^ 

Iyy(slug  sq  ft) .  790 

I  Ijj^Cslug  sg  ft) .  3  53 

LENGTH  (in)  ,  1 .  168 

CENTER  OF  GRAVITY .  100.8(0.6  1' 

distance  from  nose  (in) 

MAX.  MAJOR  AXIS  (in) .  41.7 

MAX.  MINOR  AXIS  (in) .  13.86 


Inertial  acceleration  commands  are  applied  in  polar  coordi¬ 
nates  (i.e.  ,  magnitude  of  command  applied  to  the 

pitch  autopilot  and  the  direction  (  4^0  )  applied  to  the 

roll  autopilot)  . 

The  yaw  autopilot  is  slaved  to  the  roll  autopilot  to  mini¬ 
mize  sideslip  angle  by  coordinating  the  missile  yaw  and  roll 
motion. 

Achieved  maneuver  plane  or  inertial  acceleration  in  rectan¬ 
gular  coordinates  (i.e.,  and  determined  by 

resolving  achieved  body-fixed  accelerations 

(i.e.,"!]^^  and  )  through  missile  roll  angl3(4>  )  (i.  e. , 

Euler  angles  ©•  and  are  assumed  to  be  sufficiently 

small) . 
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Figure  A- 2  Aerodynamic  Sign  Convention. 
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BTT  autopilot 


Figure  A. 3  BTT  AOIOPILOT. 
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APPENDIX  B 
LINEAR  SYSTEM  MODELS 


A.  INTHODOCTION 

The  following  appendix  addresses  the  informations  of  the 
linear  uncoupled  pitch  channel  autopilot  and  the  roll-yaw 
coupled  channel  autopilot  of  a  elliptical  airframe  taken 
from  [Ref.  2], 

The  linear  time  domain  analysis  of  the  CBTT  auropilot  was 
assumed  that  the  missile  is  initially  in  the  desired 
maneuver  plane  and  trimmed  at  ten  degrees  angles-of-attack 
(i.e.,  the  equilibrium  or  trim  angle-of-attack  in  the 

model  of  Fig.A.I  equals  10  degrees  and  the  equilibrium  roll 
rate  Pe  in  the  models  of  Fig. A,  3  equal  zero).wtien  Pe  =  0, 
Qe  (i.e.,  equilibrium  pitch  rate)  has  been  formed  to  have 
negligible  influence  in  the  lateral  model  compared 
to  and  was  therefore  set  equal  to  zero. 

B.  OHCOOPLED  PITCH  CHANNEL  AUTOPILOT 

A  general  block  diagram  of  an  uncoupled  pitch  channel 
autopilot  is  shown  in  Fig.  B.1  . 

A  normal  acceleration  command  g’s)  is  applied  to 

the  pitch  control  law  which  uses  measurements  of  missile 
body  pitch  angular  rate  (q)  and  pitch  normal  acceleration 
determine  the  required  actuator  command  (Sp^)  .  The 
actuator  is  modeled  as  a  first-order  lag  at  30Hz  (188.4 
rad/sec)  . 

The  dynamic  model  is  linearized  about  a  trim  angle-of- 
attack  as  described  in  appendix  A. 

From  the  block  diagram  of  Fig. B.1,  the  parameters  of 
interest  for  the  present  analysis  are  given  respectively 

^  C  '  Cm  *  • 
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CoNTHOL  Law  1  Actuator  Dynamic  MooHL 


Figure  B.1  Oncoupled  Pitch  Channel  autopilot. 
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For  ease  notation  the  correspondence  of  state  variables  and 
constants  are  given  in  Table  Vll 

Performing  block  diagram  manipulations  see  [aef.  5],  one 
yields  to  the  following  state  variable  equations: 

aerodynamic  equations 

XI  =  C2  A3  X2  +  C2  AU  X3  (B.1) 


X2  =  XI  -  K  Cl  ai  X2  -  K  Cl  A2  X3  (B.2) 

Control  law  equations 

X3  =  -  C3  X3  +  C3  Conv  X6  {B.3) 


X4 


Cl  C4  A1  X2  -  Cl  C4  X3  -  C4  X4 


(B.4) 


X5  =  C7  X4  -  C5  X5  -  C6  NZC  (B.5) 

Actuator  equation 

X6  =  C  E  A7  X3  “  6  X6  +  C  E  as  X8  +  C  E  A6  XII  (B.6) 


C.  COOPLED  ROLL-XiW  CHANNEL  AOTOPILOT 

A  general  block  diagram  of  a  coupled  roll- yaw  channel 
autopilot  is  shown  in  Fig. B.2. 


176 


- 

% 


CoNfROL  Law5  '-Actuators  I  Dynatiic  Model 


Figure  B.2  Coupled  Roll-law  Channel  Autopilot. 
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The  block  diagram  (Fig.B.2)  shows,  eight  parameters  of 
interest  that  are; 

-  '^ysv'.*2  =  Crp  >  «  =  A4  =  , 

A5  =  ,  »6  =  ,  A7  =  ,  and  A8  =  - 

For  ease  of  notation  the  correspondence  of  state  vari¬ 
ables  and  constants  is  given  in  Table  Vlll. 


Eqns.B.7  through  B.17  give  the  state  variable  equations 
of  the  coupled  roll-yaw  channel  autopilot: 

Aerodynamic  equations 

XI  =  C  Conv(A4  X3  +  A3  X8  +  A5  XII)  (B.7) 


X2  =  C  Conv(A7  X3  +  AS  X8  +  A6  XII) 


(B,8) 


X3  =  -  XI  -  (alpha/Conv)  X2  (B.9) 

+  KB  A  A2  X3  +  KB  A  A1  XII 

X4  =  -8  X4  -17.6  X12  ♦  17.6  PHC  (B.  10) 

Control  law  equations  (roll) 

9 

X5  =  -  (0.755/Conv)  X2  -  C  D  A7  X3  (B.  11) 

♦  (0.755  -  8  D)  X4  -  5  X5  -  C  D  A3  X8  -  C  D  A6  XII 

-  17.6  D  X12  +17.6  PHC 
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X6 


C  E  A7  X3  -  6  X6  ♦  C  E  A8  X8  «■  C  E  A6  XII 


(B.  12) 


X7  =  -  F  KC  (0.755/Conv)  X2  -  (D  ♦  E)  F  KC  A7  X  3  (B.  13) 

+  F  KC  (0.755  -  8  D)  X4  ♦  KC(15  -  5  F)  X5 

♦  KC(6F  -  15)  X6  -  15  X7  -  (D  ♦  E)  F  KC  C  A8  X3 
-  (D  +  E)  F  KC  C  A6  XII  -  F  KC  D  17.6  X12 

♦  F  KC  D  17.6  PHC 

* 

X8  =  188.4  Conv  X7  -  188.4  X8  {B.14) 

Control  law  equations  (yaw) 

X9  =  (K1  A  A2/T1)  X3  -  X9/T1  +  (K1  A  a1/Tl)  XII  (B.15) 


X10  =  (K2/Conv)  XI  -  (K2  H  ALPHAB)  /Conv  X2  (B.  16) 

+  (K2/10)  ((K1  A  A2)/T1  +  B  A4  -  H  ALPHAB  C  A7  X3 

+  (K2/10)  (B  A3  -  H  ALPHAB  C  A8  X3  ♦  K2(1  -  (1/10  T1)  X9 
+  (K2/10)  ((K1  A  A1)/T1)  ♦  B  A5  -  H  ALPHAB  C  A6)  XII 
Actuator  equations 

XII  =  188.4  conw  X10  -  188.4  XII  (B.  17) 


X12  =  X2/Coav  (B.  18) 
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D.  AEEODYHAHIC  DATA-LINEAR  APPROXIS ATION 

A  linear  approach  was  used  in  the  design  and  stability 
analysis  of  the  aultopilots  of  tha  pitch,  yaw,  and  roll 
channels,  both  uncoupled  and  coupled.  According,  a  linear 
approximation  of  the  aerodynamic  derivatives  at  m=3.95  was 
provided,  about  which  the  system  could  be  perturbed.  These 
linearized  aerodynamic  derivatives  are  given  in  [Ref-  2]  and 
are  presented  here  in  Table  IX. 
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TABLE  ?II 

Correspondence  of  Syabols  (Oncoapled  Pitch  Autopilot) 


q 


6 


p 


X 

Y 


ot 

Sp 

qs/d 

57..3qsd/I 


188.4 

150.0 

6.0 


0.53544 

0.48 

0.38375 

3.07 

0.48 

57.3 


XI 

X2 

x3 

X4 

x5 

X6 

A1 

A2 

A3 

A4 

Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

K 

Conv 
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T&BLE  VIII 

Correspondence  of  Symbols  (Couplsd  Roll-Yaw  Autopilot) 


1 


r  . 

3.  1416 

p  . 

57.3 

fi  . 

qs/d 

X  . 

qsd 1/1 2^ 

Yl  . 

qsdl/I^x 

XI  . 

0.05033 

^  ftc  . 

0.078 

. 

0.25 

Y  . 

oC  e 

he  . 

<3<^e  /57.3 

Sr  . 

0.48 

. 

0.25 

. 

. 

0.84 

6.08 

. 

1.0 

f'SR. 

. 

4.17 

C-  :  . 

0.48 

%y  . 

2.0 

. 

.  A7 

? 

. 

s 

Conv 

A 

B 

C 

D 

E 

F 

ALPHA 

ALPHAS 

KB 

T1 

K1 

K2 

KYP 

KC 

K 

D1 

QB 
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TABLE  IX 

Linearized  Aerodynamic  Derivatives 

o<  =  10° 


0.  054 

. 

c  . 

. .  -f 

0.  024  ■ 

Ca  . 

0.  027 

5(3 

c..  . 

.  ^ 

0.  015 

Y,fr 

0.  039 

. 

c«  . 

0.  010 

Cw  . 

0.  006 

. .  + 

0.  014 

. . 

. .  + 

0.  023 

Ur. 

c.  . . 

.  ^ 

0.  220 

Nc< 

c.,  . 

. . .  + 

0.  020 

0.  0137 

. 

»ne< 

c  . 

0.  055 

»Ttfp 
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aPPENDII  S 

NONLINEAR  SYSTEM  MODEL 


A.  INTRODOCTION 

A  general  block  diagram  of  the  nonlinear  model  is  shown 
in  Fig. C.1  The  three  dimensional  nonlinear  aerodynamic  model 
presented  here  is  for  the  same  conditions  used  in  the  linear 
model  shown  in  appendix  B.  The  same  flight  condition  used 
for  the  linear  case  is  the  same  for  the  nonlinear  model  (i. 
e.  ,  60  kft  altitude,  Mach  3,95).  The  control  laws  are  the 
same  used  for  the  linear  models  except  for  a  minor  modifica¬ 
tion  to  the  coordinating  branch  dependence  on  angle-of- 
attack  and  also  the  inclusion  of  anti-gravity  bias  as  stated 
in  [Ref.  2], 


B.  CONTROL  LARS 

The  control  laws  used  for  the  following  nonlinear  3-D 
studies  were  the  same  as  those  used  for  the  linear  models, 
except  for  the  gain  o<  shown  in  the  bold  line  of  coordina¬ 
tion  branch  in  Fig.  C.1.  The  new  gain  ^  is  held  constant 
at  one  degree  magnitude  for  angle-of-attack  less  than  one 
degree  positive.  For  angle-of-attack  greater  than  one 
degree  positive,  the  gain  is  equal  to  the  angle-of- 

attack.  This  keeps  coordination  for  very  small  angle-of- 
attack. 

Gravity  effects  were  not  included  in  the  linear  models 
in  appendix  B,  because  it  was  assumed  to  have  a  negligible 
influence  on  autopilot  stability  and  response  for  pertuba- 
tions  about  a  missile  trim  condition.  However,  gravity 
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Figure  C«  1 


Nonlinear  Model. 
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effects  were  included  in  the  following  nonlinear  model  where 
the  missile  body-fixed  yaw  axis  will  be  subjected  to  the 
full  force  of  gravity  and  may,  therefore,  have  a  significant 
influence  on  sideslip. 

The  gravity  effects  are  compensated  for  by  pitch  and 
roll  motion  of  th^  missile  which  have  less  influence  on 
sideslip  than  yaw  motion.  In  inertial  rectangilar  coordi¬ 
nates  one  has: 

,  ^  =acceleration  command  in  inertial  direction 


=  V|^^_  cos(e) 


(C.1) 


where 

*'12.  =  guidance  command (  g’s  ) 

-Cos©  =  anti-gravity  bias  command  (  g's  ) 

=  acceleration  command  in  the  inertial  direction 

guidance  command (  g's  ) 

There  is  no  gravity  effect  in  Y  direction.  In  order  to 
assure  an  anti-gravity  bias  of  just  one  g,  oa e  needs  to 
modify  the  anti-gravity  command  as  follows; 

elliptical  airframe  inertial  acceleration  command 

=  Nzc-0. 91  3cos  (^) 

c 


C.  PEBFOHANCE  COHHANDS 

The  commands  2  gees  is  first  applied  in  the  0°  or 
upward  direction  at  2  seconds.  since  both  the  missile 
roll  angle  and  roll  angle  command  are  at  zero  degrees,  there 
is  no  roll  motion  and  the  missile  turns  upward  as  a 
skid-to-turn  controlled  missile.  At  5  seconds,  a  second 
2  g's  inertial  guidance  command  is  applied  in  the  downward 
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or  1800  direction.  The  missile  is  commanded  to  roll  through 
1800  while  mowing  in  a  coordinated  manner  in  yaw  and  roll  to 
minimize  sideslip  angle  and  prevent  or  minimize  negative 
angle -of -attack. 

The  response  of  achieved  maneuver  and  cross-plane  accel¬ 
erations  during  the  first  2  seconds  are  due  to  initial 
conditions,  gravity  and  anti-gravity  bias  effects.  The 
initial  conditions  were  added  to  minimize  the  transients 
which  result  when  gravity  bias  commands  the  autopilot  for 
constant  altitude  missile  flight.  The  initial  conditions 
were : 

^  =  angle-of-attack  =  2.41  deg 

Sjp  =  pitch  tail  angle  =  0.658  deg 
©  =  pitch  Euler  angle  =  3.65  deg 

output  of  pitch  acceleration  feedback  lag  =  -1.0  g*s 
pitch  control  law  acceleration  error 
lag  prior  to  dc  gain  =  -  0.0105 
5pc  =  pitch  actuator  command  =  0.558  deg 

The  achieved  maneuver  plane  acceleration  (Eqn.C.2)  is 
calculated  from  the  body-fixed  acceleration 
the  roll  angle  {)  as  follows; 


^6 


and 


"I 


(C.2) 


During  the  first  command,  achieved  body-fixed  yaw  accel¬ 
eration  and  missile  roll  angle  are  equal  to  zero, 

because  the  roll  channel  is  not  commanded.  Therefore, 
achieved  maneuver  plane  acceleration  is  equal  to  the  body- 
fixed  acceleration  y[  z^.  During  the  second  command  the 
missile  roll  angle  has  the  same  roll  angle  response  of  the 
roll  autopilot. 
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D.  AERODYHAMIC  DATA  -  NONLINEAR  REPRESENTATION 

The  aerodynamic  data  are  given  in  [Ref.  2].  The  entire 
study  was  conducted  using  (!=3.95  and  the  aerodynamic  coeffi¬ 
cients  are  based  on  a  body- fixed  axis  system  of  Pig.B. 2. 

For  reference,  the  normal  force  and  pitching  moment 
plots  given  in  appendix  B  of  [Ref.  2 ]  at  8=3.  95  were  repro¬ 
duced  in  Table  X  and  XI,  which  give  just  the  values  for 
points  of  interest.  These  tables  present,  too,  the  corre¬ 
spondent  derivatives  that  were  found  by  simulation  using  a 
fortran  program  shown  in  appendix  F.  The  aerodyiamic  deriv¬ 
ative  ofCjj,  ,  and  Cy  with  respect  to  sideslip 

angle  (h  ,  yaw  control  6y  ,  and  roll  control  are 

presented  in  Table  XII,  XIII  and  XI7  which  give  just  values 
for  points  of  interest  and  includes  the  correspondent  deriv¬ 
atives  that  were  found  by  simulation  using  a  fortran  program 
given  in  appendix  P. 

The  block  diagram  of  Fig.C.1  shows  10  parameter  of 

interest  that  are  : 

A1  =  C^(oc,<fp)  ,  A2  =  ,  A3  =  , 

A4  =  C  (X)  ,  A5  =  C.  (X)  ,  A6  =  Cy  (x)  , 

A7  =  »  A8  =  (X)  ,  A9  =  , 

and  AA  =  Cy^  (<^)  . 

Using  the  same  procedure  as  given  in  appendix  B  one  can 
have  the  Tables  X7  and  XVI  that  give  the  correspondence  of 
symbols  for  the  nonlinear  model.  Egns.C.3  through  C.21  give 
the  state  variable  equations  of  the  nonlinear  model. 
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XI  =  -  C4  XI  -  C02  C4  CN(X6,X4)  (C.3) 


X2  =  C7  XI  -  C5  X2  ♦  C6  C19  C0S(X17)  -  C6  NZC  (C.4) 


X3  =  -  C7  C8  X1  ♦  (C5  C8  -  C9)  X2  +  (C9/Conv)  X5  (C.5) 


+  C6  C8  NZC 

# 

X4  =  -  C3  X4  ♦  C3  Conv  X3  (C,6) 

« 

X5  =  (X17  Xl8)/Conv  ♦  Cl  Conv  CM(X6,X4)  (C.7) 

X6  =  X5  -  KB  CM(X6,X4)  -  (X16  X18|/Con7  (C.8) 

♦ 

X7  =  X5  Cos(Xl9) /Conv  -  X17  Sin  (XI 9 ) /Conv  (C.9) 

f 

X8  =  -  C10  X8  +  (C12  -  C11  C13)  X10  (C.  10) 


XI  5 


(C.  11) 


♦  C11  C14  PHC  -  C11  C14  X19  -  CO  C1l(Cy  (X6) 

♦  Co^  (x6)  X13)-  C12  X18/Conv  ' 

X9  =  -  X9/T1  +  (K1/T1)  C02(C  (X6)  X16 


♦  C^^(X6)  XI 5) 

T^y 
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X10 


C13  X10  ♦  C14  PHC 


C14  X19 


(C.  12) 


XII  =  -  C15  XII  ♦  C16  C0(C^^(x6)  X16  (C.13) 

♦  C  (X6)  X15  ♦  C  (X6)  X13) 

X12  =  -  C7  X12  ♦  K  (C17  -  C10/C18)  X18  (C,  14) 


+  (K/C18)  (C12  -  C11  C13)  X10  +  K  C11  (C14/C18|  PHC 

-  K  C11  C14/C18  X19  -  (K/C18)  CO  C  (X6)(C11 

+  C16)  X16  -  (K/C18)  CO  Cjp(X6)(C11  +  C16)  X15 

-  (K/C18)  CO  Cj{^(X6)  (C11  +  C16)  X13 

-  (K  C12/CDnv)  XI 8  +  {C15/C18  -  C17)  XII 

# 

X13  =  -  C3  X13  ♦  C3  Conv  X12  (C.  15) 


X14  =  (K2/10)  (-  X9/T1  +  K1  T1  C02  (C  (X6)  X16  (C.  16) 

+  Cy^^(X6)  X15)  -  (X15  XI  8) /Conv2  (X5  -  KB  C02  CN(X6rX4) 

-  (X16  X18)/Conv  X18  -  (KXP/Conv)  (X5  -  KB  C02  CM(X6,X4) 

-  (X16  X18)/Conv  X18  -  (KYP/Conv)  CO  XX6  (C^^(X6)  X16 

-  Cjj^{X6)  X15  *  C,^p(X6)  X13)  ♦  (K2/Conv)  X17 

+  K2  X9  -  (K2  KYP/Conv2)  XX6  X18 


X15  =  -  C3  X15  ♦  C3  Conv  X14 


(C.  17) 


X16  =  KB 


C02  (Cy^  (X6) 


X16  * 


(X6) 


X15 


+  (X16  X18)  /Conv  -  X17 


(C.  18) 
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X17  =  -  (X5  XI 8) /Conv  ♦  C01  Conv(3^^(X6)  X16 
+  (X6)  X13) 

X18  =  CO  Conv  (Cn  (X6)  X16  +  C^  (X6»  X15 

*  C  (X6)  X13) 


XI 9  =  X18/CDnv 


(C.  19) 


(C.  20) 


(C.  21) 
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T&BLE  X 

Normal  Force  Coefficient  and  Derivatives  I!=3.95 


^  Cu 

^  C  )s4  <fp) 

1)  <£|^ 

-4.00C0 

-lO.OOOC 

-1.C500 

0.1842 

0.0308 

-4. COCO 

0.0000 

-0.8000 

0  .1453 

0 .0208 

-4.C0CC 

5.000C 

-0.7000 

0.2496 

0.0196 

-4.0000 

10.0000 

-0.6000 

0.2141 

0.0208 

0.0000 

-10.0000 

-0.2500 

0  .2110 

-0.0332 

0.0300 

o.oooc 

-o.caoo 

0.20  74 

0.0448 

0.0000 

5.0000 

0.1300 

0.1827 

0.0336 

0.0000 

lO.OOOC 

0.2000 

0.1929 

-0.0112 

4.0000 

-10.0000 

0.6000 

0.2092 

0.0200 

4.0000 

0.0000 

0.8000 

3 .2253 

0.0200 

4. COCO 

5.000C 

c.sooo 

0.2196 

0.0200 

4.0000 

lO.OOOC 

l.COOO 

0.2141 

0.0200 

8. COCO 

-lO.OOOC 

1  .4500 

0  .2272 

0 .0042 

8.0000 

O.OOOC 

1.7000 

0.2265 

0.0342 

8.0000 

5.0000 

1.6500 

0.2288 

0.0229 

8.0000 

lO.OOOC 

1.9000 

0.2256 

-0.0058 

12.0000 

-lO.OOOC 

2.4000 

0.2321 

0.0200 

12.0000 

0.0000 

2.6000 

0  .2188 

0.0200 

12. COCO 

5.000C 

2.7000 

0.2153 

0.0200 

12.0000 

lO.OOOC 

2.6000 

0.2335 

0.0200 

16.C0CC 

-lO.OOOC 

3 .3000 

0 .2319 

0  .0083 

16.0000 

0. OOOC 

3.50C0 

0.2482 

0.0283 

16.0000 

5.0000 

3.6500 

0.2599 

0.0308 

16.0000 

lO.OOOC 

3.6000 

0.2654 

0.0283 

20.0000 

-10.0000 

4.3000 

0.2653 

0.0300 

20. COCO 

o.oooc 

4.6000 

0.2884 

0.0300 

20.0000 

5.0000 

4.7500 

0.2825 

0.0300 

20.0000 

10.0000 

4.9000 

0 .2798 

0.0300 

24.0000 

-lO.OOOC 

5.40CO 

0.2819 

0.0300 

24.0000 

0.0000 

5.7000 

0.2482 

0.0300 

24.00CC 

5.0000 

5.6500 

0  .2599 

0.0300 

24.0000 

lO.OOOC 

6.C000 

0.2654 

0.0300 
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■V 


TiBLB  XI 

Pitching  Moment  Coefficient  and  Derivatives  (!=3.95 


o< 

S 

ip) 

'Doc 

-4.000C 

-lO.OOOC 

0.4850 

0.0130 

-0.0159 

-4.0000 

-5.0000 

0.2800 

-0.0200 

-0.C61C 

-4.0000 

0.0000 

-0.C750 

0 .0190 

-0.0759 

-4.00C0 

lO.OOOC 

-0.^300 

0.C197 

-0.0149 

0.0000 

-10.0000 

0.5500 

0.0185 

-0.055C 

O.OOCC 

-5.000C 

0.2750 

0 .0131 

-0.0550 

0.0000 

0. 0000 

-o.cooo 

0.0186 

-0. 0550 

0.0000 

10.0000 

-0.5500 

0 .0195 

-0.0550 

4.00CC 

-lO.OOOC 

0.6250 

0.0180 

-0.0551 

4.0000 

-5.000C 

0.2500 

0.0200 

-0.055C 

4.0000 

0.0000 

0.C750 

0 .0190 

-0.0551 

4.0000 

lO.OOOC 

-0.4800 

0.0147 

-0.0561 

8.000C 

-10.0000 

0.7100 

0.0294 

-0.0563 

8. 00  CO 

-5.000C 

0.^300 

0.0232 

-0.0558 

8.0000 

O.OOOC 

0. 1500 

0.0178 

-0.0563 

8.0000 

10.0000 

-0.4300 

0.0118 

-0.0603 

12.000C 

-lO.OOOC 

0.6500 

0.0331 

-0.0640 

12.0000 

-5.0000 

0.5250 

0.0184 

-0.0655 

12.0000 

0.0000 

0.2000 

0  .0034 

-0.0640 

12.0000 

10.0000 

-0.3900 

0.0057 

-0.052C 

16.0000 

-10.0000 

0.9500 

0 .0182 

-0.0725 

16.0000 

-5.000C 

0.5750 

0.0120 

-0.0763 

16.0000 

0.0000 

0.2000 

0.0059 

-0.0725 

16. COCO 

lO.OOOC 

-0.4000 

-0  .0122 

-0.0425 

20.0000 

-lO.OOOC 

l.CCOO 

0.0065 

-0.0756 

20.0000 

-5.0000 

0.6220 

0 .0065 

-0.0755 

20.00CC 

O.OOOC 

0.2450 

0.0067 

-0.0752 

20.0000 

10.0000 

-0. 5000 

-0.0395 

-0.0736 

24.00C0 

-10.0000 

1 .COOO 

-0  .006  8 

-0.0841 

24.0000 

-5.000C 

0.5870 

-0.0290 

-0.0818 

24.0000 

0.0000 

0.1750 

-0.0516 

-0.0837 

24.00CC 

lO.OOOC 

-0.7300 

-0.0772 

-0.1001 
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TABLE  XII 

Sideslip  Derivatives  and  Derivatives 


-4.00C0 
0.0000 
4.0000 
5.20CC 
8.00C0 
10.0000 
12.60CC 
16.0000 
20. COCO 


0< 


-4.00CC 
0.00  00 
4.0000 
5.20CC 
8.0000 
10.0000 
12.6000 
16.0000 
20.00CO 


0.0240 
0. C240 
0. 0240 
0.0240 
0. C240 
0.C240 
0. C270 
0. C290 
0.0320 


Cy  {.e^  ) 


•0. C490 
•0. 0440 
-0.0480 
•0. C490 
•0.0520 
•0. C550 
■0. C580 
•0. 0610 
•0. C650 


«a< 

Cn('-) 

V 

""t)  oc 

-4.0000 

0.0110 

-0.0027 

0.0000 

0.0000 

-0.00  27 

4.00CC 

-0.  Clio 

-C.  0022 

5.2000 

-0.0135 

-0.0024 

8. COCO 

-0. 0220 

-0.0027 

10.0000 

-0.  0270 

-0.0021 

12.6000 

-0.0310 

-0.0015 

16.0000 

-0. 0360 

-0.0  014 

20.0000 

-0.0410 

-0.0011 

C< 


C.OOOO 
C.0000 
G.OOCO 
0.0000 
C.00C0 
C.0005 
0.0  009 
C.  0  0C7 
0.0008 


^  Cy 


c>4 


0.0  024 
0.0001 
-0.0009 
-C.OOC9 
-0.0013 
-0.0013 
-C.0010 
-0.0009 
-0.0011 
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TABLE  XIII 

Yaw  Control  Derivatives  aod  Derivatives 


oC 


-4.00CC 
0.0000 
4. 0000 
5.20C0 
8.0000 
10. COCO 
12.60CC 
16.0000 
20.0000 


0. 0150 
0.0150 
0.0150 
0. 0150 
0.0150 
0. 0150 
0. 0150 
0.0170 
0. 0195 


D  wL 

-4.0000 

0. 0025 

-C.  00C7 

0.0000 

0.0000 

-0.0006 

4.00CC 

-0. C020 

-0.0  014 

5.2000 

-0.0040 

-0.0  015 

8.00C0 

-0.0075 

-0.0013 

10. COCO 

-0.  ClOO 

-C.0  014 

12.6000 

-0.0140 

-0.0015 

16.C000 

-0. C190 

-0.0014 

20.0000 

-0. C240 

-0.0011 

o< 

Ch 

^  oc 

-4.00C0 

-0. 0400 

0.0000 

0.0000 

-0. 0400 

C.0  0C0 

4.0000 

-0.0400 

C.OOCO 

5.20CC 

-0.  C400 

0.0000 

8.0000 

-0.0400 

C.00C0 

lO.OOCO 

-0. 0400 

C.OOCO 

12.600C 

-0. C400 

-C.00C3 

16.0000 

-0.0420 

-0.0007 

20.000C 

-0.  C453 

-0.0010 

u. 


0.0000 
0.0000 
0.0000 
0.0  000 
0.0000 
0.0000 
0.0003 
C.0006 
0.0  006 
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TABLE  XIV 

Roll  Control  Derivatives  and  Derivatives 


o< 


-4.00CC 

O.QOCC 

4.0000 

5.20CC 

8.00CO 

10.0000 

i2.eocc 

16.0000 

20.0000 


op. 


•0.0045 

0.  cooo 

0.0045 
0. C055 

0. Clio 

0.0145 
0. C195 
0.0250 
0. C320 


-4.0000 

0.0232 

C.00C0 

O.COCO 

0. 0232 

0.0  0  00 

4.0000 

0. 0232 

C.00C0 

5. 2300 

0. 0232 

C.  0  000 

8.0000 

0.0232 

0.0  0  CO 

10.0000 

0.  0232 

0.0  003 

12.6000 

0. 0250 

0.0  0  06 

16.0000 

0. 0270 

0.0006 

20.0000 

0.  C298 

C.  0  0C8 

"~d  oc 


Q.OOll 
0.0  011 
C.00C9 
0.0012 
C.  0  018 
0.0  018 
C.  0  018 
0.0017 
0.0018 
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TABLE  XV 

Correspondence  of  Symbols  (  nonlinear  Model  ) 

. XI  p  . 

T 

.  X2 

<p  . 

.  X3 

. 

. X4 

. 

q 

. X5 

. 

of. 

.  X6 

B 

. X7 

1  Sy 

. . 

T1 

.  X8 

. 

Y2 

. X9 

'I?;  . 

XI 

.  X10 

>  . 

^  . 

X2 

. XII 

Si^c 

. X12 

. 

Sr. 

. X13 

qSd/Ip^^  •••••••«• 

. X14 

qSd/Iz2  . 

S-y 

. X15 

qSd/W  . 

P> 

.  X16 

r 

. X17 

3.1416  . 

J 
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TABLE  XVI 

Correspondence  of  Symbols  (Coat.  Nonlinear  Sodel) 


188.4  . 

150.0  . 

6.0  . 

0.530544  . 

0.48  . 

. C7 

0.38375  . 

3.07  . 

5.0  . 

0.05033  . 

0.755  . 

8.0  . 

17.6  . 

6.0  . 

0.078  . 

15.0  . 

4.0  . 

0.913  . 

0.25  . 

0.839  . 

6.08  . 

1.0  . 

'~c> 

'li  ot 

^  c< 

C|^«  Cfjt/ 


DA16 


DA14 


DA26 


DA2  4 


DA3 


DA4 


DAS 


DA6 


DA7 


DAB 


DA9 


DAA 
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UNCOUPLED  PITCH  AUTOPILOT  -  CSMP  PROGRAM 
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AERODYNAMIC  EQUATIONS 
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COUPLED  ROLL^YAW  AUTOPILOT  -  CSMP  PROGRAM 
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